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by the Space Systems Division of LMSC, under the management of the NASA Project 
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Section 1 
SUMMARY 


A thermodynamic vent system is one wherein vent fluid is extracted from a tank, 
expanded through a Joule -Thompson valve to produce a temperature drop, and then 
passed through a heat exchanger to cool the remaining fluid in the tank. The program 
described in this report was designed to evaluate a passive system for which the heat 
exchanger is formed by wrapping tubing over the tank wall. No active mixing is used 
to stir the tank contents. Therefore, conduction and/or free convection will control 
the energy distribution within the tank and the rate at which the tank pressure can be 
controlled by cooling at the tank wall. This suggests that the system efficiency, 
which is based on the ability to control pressure while minimizing total mass loss, 
may be gravity dependent. 

The objectives of this program were to: (1) develop analytical models which describe 
the mass and energy distribution inside cryogenic tanks and relate these to the pres- 
sure control characteristics (primary variables were gravity level, heat flux,- and 
tank size), (2) develop design models for the wall mounted system, (3) devise and 
implement a one-g test program to simulate the low-g conditions and provide data 
to supplement the analytical models, and (4) evaluate the validity of the models. 

The program was divided into the following five technical tasks: 

Task I — . Characterization of Cryogenic Tank Pressure Histories 

Task II — Experimental Evaluation of the Heat Transfer Model 

Task HI — Design and Fabrication of System for Experimental Evaluation 
of Design Model 

Task IV - Test Evaluation of Design Model System 
Task V - Correlation of Analytical and Experimental Data 
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Two models were developed in Task I; one for the case of pure conduction (zero-g) 
and one for the case where heat transfer was convection dominated. Further, the 
zero-g case requires that a model be developed for each tank shape (i.e. , sphere, 
cylinder with L/D = 2, 3, 4, etc.). An R-C network is generated for each con- 
figuration and ullage volume with the liquid -vapor interface assumed to be spherical. 
Phase change occurs anywhere in the liquid if the local temperature is above the 
instantaneous saturation temperature. The R-C nodal network is put into a thermal 
analyzer program along with criteria for phase change. The program computes the 
mass and energy balances and corresponding pressure and temperature histories 
for a given heat exchanger resistance (wall-to-tube) . This model indicates that the 
tank wall heat exchanger will control tank pressure, even under zero-g conditions. 

With an active mixer system the tank contents are maintained in a funiform condition 
with mechanical mixing. Thus, there is no subcooling in the bottom of tank which 
results in excess loss of vented propellant. With the tank wall heat exchanger colder 
fluid collects in the bottom during the cooling cycle giving rise to nonuniform conditions; 
this results in some subcooling and a concommitant excess vented propellant. However, 
even in zero-g this becomes a relatively insignificant penalty when the vehicle has 
been through 5 or 6 vent cycles. 
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Section 2 
INTRODUCTION 


So effective have become the thermal protection systems for cryogenic space pro- 
pulsion vehicles that nonvented storage of liquid cryogens is now a practical considera- 
tion for Earth-orbital and planetary orbit injection missions. Even so, in all pro- 
pellant tankage systems a pressure relief system is mandatory for safety and system 
operational considerations. For example, venting maybe desired for reduction of 
tank pressure after an engine firing. In this way, lower tank operating pressures 
and tank weights may be used and still provide the required net positive suction 
pressure to the engine, when successive firings are demanded by the mission profile. 

Normally, a gas-vent relief valve might be considered for venting the cryogen tank. 

In space flight, however, long periods of near zero-gravity coast of the spacecraft 
occur. Space flight experience has shown that in zero-gravity the location and move- 
ment of the liquid propellant in the tank is uncertain. Thus, under such conditions, 
pressure relief of the tank through an ordinary relief valve is unreliable and possibly 
very inefficient because of the likely ingestion of large amounts of liquid cryogen that 
would be vented directly overboard. Actual mission failures have occurred beoause of 
excessive propellant loss and consumption of attitude control propellants to correct 
for large variations and unbalances in vehicle motions that were induced by the vent 
system. 

One concept for providing effective control without the above mentioned problems is 
a thermodynamic vent system. During operation of this system, vent fluid is with- 
drawn from the tank and expanded through a Joule -Thompson valve to produce a 
refrigerated vapor. This cold fluid is then passed through a heat exchanger that is 
in contact with the tank propellants, where it absorbs enough energy to vaporize the 
liquid in the wet vapor. The vented propellant is then vented overboard as a saturated 
or superheated vapor just as though it had been extracted directly from the ullage. 

The thermodynamic process is illustrated on figure 2-1. 
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Fig, 2-1 Thermodynamic Process 

The NASA-LeRC recognized the potential of this type of system some time ago and 
awarded contracts to LMSC (NAS 3-7942, 3-12033) which resulted in the design, 
development and experimental demonstration of an active system that used a compact heat 
exchanger. 


In this previous work, design data and analysis techniques were developed which 
related both the thermal and hydrodynamic (mixing) characteristics of the vent system 
to the pressure response in the tank. It was shown that when the mixer jet velocity 
exceeded a critical valve necessary to give complete circulation, the tank pressure 
could be predicted with mixed model theory and pressure response was limited by the 
thermal capacity of the heat exchanger. The analytical models were confirmed with 
approximately 1500 hours of testing in 41.5 inch and 109 inch diameter tanks. The 
results of this previous work on the active system are presented in Ref. 1-3. A similar 
program is described in Ref. 4. 

The basic objective of this program was to provide techniques and data for the passive 
tank wall heat exchanger concept that would parallel that previously obtained for the 
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active system. The program is both analytical and experimental in nature and is 
designed to measure the effectiveness of the tank wall heat exchanger system as a 
low gravity tank pressure control device. Specific objectives which served to accom- 
plish the foregoing broad objective were as follows: 

1. Develop analytical models which describe the energy distribution inside 
cryogenic tanks and which relate these to the pressure control characteristics. 

2. Develop design models for the wall mounted heat exchanger. 

3. Devise and implement a one-g test program for simulating low-g convection 
and provide data to supplement the analytical models. 

4. Correlate the experimental results with the analytical models. 

Performance curves have been generated for a wide range of tank sizes, g-level, heat 
flux, etc., and are contained herein. However, the theoretical developments, evaluation 
criteria, and test techniques can be used for any specific application, such as a cryogenic 
space tug. They are also valid for non-cryogenie fluids. Details of the program 
are described in the following sections of the report. 
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Section 3 


CHARACTERIZATION OF CRYOGENIC TANK PRESSURES 


When no active mixing is used to stir the contents, conduction and/or free convection 
will control the energy distribution within a space vehicle cryogen tank. The extent 
to which one mode will limit the pressure control process depends upon the heat flux 
and g-level. This is illustrated on Fig. 3-la as a thermal resistance dependency. If 
the g-level is sufficiently high to induce free convection it takes time to establish a 
boundary layer. During this transient period, the thermal resistance could be 
expected to vary in a manner depicted in Fig. 3 -lb. An analysis was performed to 
estimate the time required to develop the boundary layer. However, the approach 
in this program was to determine the pressure response using conduction models, and 
convection model with fully established flow. If the convection model results in vent 
periods of the same order, or shorter, than the predicted time for developing the 
boundary layer, then the more conservative conduction model would be used for the 
heat exchanger design. 


The basic configurations and operational parameters for this study are listed in 
Table 3-1. However, this wide matrix can more easily be assimilated reducing it 
to physically meaningful dimensionless parameters, such as the Bond number (Bo) 
and Raleigh number (Ra). These dimensionless parameters indicate the importance 
of liquid orientation, free convection, conduction, etc. and thus the type of analytical 
models that are needed. 
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Thermal Resistance 



(Heat F3 ux)x(g-Level) 

(a) Heat Transfer Domains 



Time 


(b) Thermal Resistance Related to Time 

Fig. 3-1 Pictorial Illustration of Heat Transfer Domains and Thermal Resistance 
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Diameter 
m (feet) 

1.22, 2.44, 12.2 
(4), (8), (40) 

Cylinder 0.90 L/D - 2 

(2.95) 

1.80 L/D = 2 
(5.9) 

1.37 L/D = 4 

(4.5) 

3.66 L/D = 4 

( 12 . 0 ) 

8.84 L/D - 2 
(29.0) 

6.86 L/D = 4 

(22.5) 


Configuration 

Sphere 



Acceleration 

g /g Q 

0 , 10 " 6 , 1 


Heat Flux, 
w/m 2 

(Btu/Hr-Ft 2 ) 


0.31, 1.56, 3.14 
(0.1), (0.5), (1.0) 

0.31, 1.56, 3.14 


0.05, 0.50, 0.90 0, 10" 6 , 1 0.31, 1.56, 3.14 

0.05, 0.50, 0.90 0, 10“ 6 , 1 0.31, 1.56, 3.14 

0.05, 0.50, 0.90 0, 10 -6 , 1 0.31, 1.56, 3.14 

0.05, 0.50, 0.90 0, 10” 6 , 1 0.31, 1.56, 3.14 

0.05, 0.50, 0.90 0, 10~ 6 , 1 0.31, 1.56, 3.14 















The liquid -ullage interface configuration can be inferred from the Bond number 
which can be computed from the following relationship: 


Bo = aD 


2 


4<rg, 


where: a. = g/g Q » P - liquid density, a = surface tension 


g Q = standard gravitational acceleration 


g = universal gravitational constant 
c 


This Equation is graphed as Figs. 3-2 and 3-3 for hydrogen and oxygen, respectively. 

Also shown are the points from Table 3-1, identified by tank number, and g level 

(except where g = Bo = 0). As indicated on the figures, at 10 g's, the four largest 

tanks will have interface configurations somewhere between spherical and flat (i.e., 

2 

1 s Bo ^10 ). The other tanks and g-levels can be treated with either a flat or 
spherical interface model. 


The heat transfer domains to be modeled can be characterized by a modified Raleigh 
number defined by: 


* 

Ra 



This Equation is graphed on Figs. 3-4 and 3-5 for LH 0 and LO qJ respectively, which 

_ g 

also show the computed points from the Table 3-1 matrix for a = 1, and 10 . A 

$ 

non-dimensional heat flux, q^ = q^/ 1.0 has been used for convenience in plotting 
Figs. 3-4 and 3-5. The characteristic length used in the computation was approximated 
from the tank length (L) and ullage ratio (V u ) as X = L (1 - V ). 
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The study points plotted on Figs. 3-4 and 3-5 can be compared with the types of 

* 

convection indicated for various values of the Raleigh number (Ra), as shown along 

the right hand edge of these graphs. The study points lie within rather narrow bands 

~6 

along each of the acceleration levels considered, a - 0, 10 and 1, since a limited 
range of q w was considered, O.lsq^sl.O. For the one gravity study points, 

10 _1 ^ aq ?1, the flow regime should vary from transition flow between laminar 
and turbulent boundary layers for the smallest scale tanks up to fully established 
turbulent convection for all larger scale tanks. Even the smallest scale tanks with 
one gravity acceleration and a ground hold heat rate would all experience turbulent 
free convection boundary layers . 

The study points for at- 10 acceleration level and 0.1 ^ q^ * £ 1*0 wall heat flux, 

range from gross viscous motion for the smallest dimensions X to a turbulent free 

convection boundary layer for large scale X dimensions. Notice that the first five 

sets of study points, designated by tank numbers 1, 4, 2, 5 and 7, at the smallest 

X lengths, all apply to low Bond numbers, B q < 1.0, and low fill levels, (1 - *) - 0.1. 

The liquid configuration in these tanks will result in longer values of X than those es ti- 
le 

mated in this study, so that the values of Ra might be one or two orders of magnitude 
higher than shown for these points. Hence, all of the ct = 10 acceleration level study 
points would be in the laminar and turbulent boundary layer regimes . Based upon the 
Figs. 3-4 and 3-5 study points, therefore, gross viscous motion was not modeled. 

* 

For zero gravity, (a = 0), the modified Raleigh number Ra = 0. Hence, one third 
of the total study points, are in the pure conduction heat transfer domain. Although 
this acceleration level could never be practically obtained, its heat transfer model 
represents a limiting solution of the tank pressurization time. 
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TANK DIAMETER - feet 


Fig. 3-2 Bond Number as a Function of Acceleration and Tank Size for LH 

3-6 

LOCKHEED MISSILES 8t SPACE COMPANY 





BOND NUMBER - dimensionless 



TANK DIAMETER - feet 


Fig. 3-3 


Bond Number as a Function of Acceleration and Tank Size for LO 


2 
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CHARACTERISTIC LENGTH - feet 


Fig. 3 


4 Raleigh Number as a function of Length and Acceleration for LH, 
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MODIFIED RAYLEIGH NUMBER - dimensionles 




3.1 TRANSIENT BOUNDARY LAYER ANALYSIS 


One of the important considerations involved in selecting a heat transfer model for 
low-g conditions is the amount of time required to establish a natural convection 
boundary layer. Expressions have been developed by Schwartz and Adelberg (Ref. 5) 
for such time constants for the case of uniform wall temperature and by Siegal (Ref, 6) 
for a laminar boundary layer with uniform heat flux. An analysis has been conducted 
to derive the time constant for the turbulent boundary layer with uniform heat flux. 
The system used in the analysis is illustrated in Fig. 3=6. 



Fig. 3-6 Model for Boundary Layer Analysis During Transient 
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Vliet (Ref. 7) gives the following expressions for boundary layer thickness and for 
a fictitious velocity which characterizes a fully developed, turbulent boundary layer 
with uniform heat flux. 


U x = 13.15 v 



5/14 


5-14 


,'Vt 


3/7 



1/14 



1/14 


5/17 

x 


Using the velocity profile from Ref. 8, the steady state velocity at any distance y 
from the wall, in the X direction is given by: 


U s /U 1 = 1/7 |l - y/sj 4 

The average velocity in the boundary layer is obtained by integrating over the 
thickness of the boundary layer. Therefore, 

U = i/fi f U dy = 0.147 U 

S S J 5 

o 


By substitution, we get: 


U 

s 


= 1.93 v 




3/7 

x 


Adelberg (Ref. 5) used a Taylor series to describe the thickness ratio between the 
transient and fully developed boundary layer, i.e., 
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g (x,t) 

fi s 


= 2 ff) ■ ff f f° s * 

vs/ V s/ \ 



where: 6, = thickness at location x, and time 

(x,t) 

t = time 

t = time required to reach steady state 


The equations for 6 g 


and TT 

s 


can be combined to give: 


U 

s 



Assuming the expression for U , to hold at any time, we can get by integration: 

s 



Conservation of mass in a control volume of unit width gives: 


t t 

s s 

/(M e -M o ) dt =/ M STORED dt = P5 s L 
o o 


where: 


M — PU, 6 
o (x,t) (x, t) 


o 1 



(Ref. 1) 


and t is a reference time constant which must be significantly less than t . 
^ s 
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Substituting the thickness ratio into the expression for U , , we get: 

A 9 L 


t t 

|M o dt =/ M s [ 2 (^s) - (*A.) 


8/5 


dt =0.57 M t 
s s 


o 


Combining the two preceding equations, gives: 




0.43 M g t g = P6 g L . 1 


At steady state conditions the mass flow rate M is proportional to the average 

s 

velocity and the boundary layer thickness at any location x . Evaluating at x = L, 
we have 


M s =P& S> L 6 S,L 


Combining the preceding two expressions, gives: 


°.4SP U S)L <5 SiL t s =M s L 


T / 6 s dX = ? / 12 fi s,L 


These equations combine to give an expression for t in terms of the fully developed 

S 

average exit velocity and the boundary layer run length. 


t g = 1.37 


( L/C s.l) 
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The bracketed quantity is the time it would take a particle of fluid moving at the 
fully developed velocity U T to traverse the length of the boundary layer. During 

S y ±J 

the transient, this time is effectively increased by 37 percent. 


t 


s 


0.71 

v 


(pv) 5/M 


J 



5/14 


Figures 3-7 and 3-8 show the calculated transient for hydrogen and oxygen, 
respectively. Also shown is Sie gal’s solution for laminar boundary layers. Sie gal’s 
solution gives: 


t 


s 


4 . 33 (1. 68 + Pp) 


2/5 



-2/5 p 2/5 x 2/5 

r 

a 


The resulting transient times can be compared to the vent cycle period as determined 
with a pure conduction model and with a steady state, free convection model to gain 
insight into the controlling mechanism. It can be seen on Figs. 3-7 and 3-8 that for 
the large tanks, the time constants are on the order of 10-20 hours. 


3.2 ZERO-GRAVITY MODEL DESCRIPTION 


The zero-gravity model of the propellant tanks and thier contents is based upon 
conduction heat transfer between the tank, ullage, and liquid nodes, with vaporization 
or condensation changes of phase of the propellant. The liquid was considered to 
have a fixed zero-gravity orientation in the tank, and the liquid volume was divided into 
nodes. These liquid noses and the tanks wall nodes were then programmed into THERM, 
which computes the heat transfer in the tank system. 

Once a nodal structure has been generated for a given geometry (L/D) any tank 
diameter can be run by inputing the appropriate scaling factors, and any fluid can 
be used by inputing the appropriate fluid properties. 
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Fig. 3-7 Transient Periods for Development of Boundary Layers in Liquid Hydrogen (q^ = constant) 
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Transient Periods for Development of Boundary Layers in Liquid Oxygen (a = constant) 



Liquid/ullage coupling model was programmed in subroutine external to THERM, and 
at each THERM time step, new values of ullage volume, ullage vapor mass, ullage 
temperature, and tanks saturation temperature and pressure are computed based upon 
the THERM mode temperatures and heat rates. 

3. 2. 1 ULLAGE Subroutine 

This subroutine is used to solve a matrix of four simultaneous equations which relate 
energy and mass balances of the liquid and vapor in the propellant tank. 

Primary effort in this model was spent in simulating the evaporation and condensation 
effects that would occur in zero-gravity conditions. Besides evaporation or condensa- 
tion that would always occur on the normal liquid/vapor interface; vaporization and 
subsequent condensation is allowed on the liquid covered tank nodes or in the bulk liquid 
nodes, while condensation and subsequent evaporation is allowed on the normally dry 
tank wall nodes opposite the ullage space. Liquid or vapor respectively formed by 
condensation or evaporation is always assumed to remain at the same location where 
it was formed, and the mass of such liquid or vapor was continually summed at each 
tanks wall node and bulk liquid node. The opposite and subsequent change of phase is 
then limited by the mass formed originally, and if this mass was exhausted, the 
corresponding node is then returned to its normal thermal capacity status in THERM. 

In this way, the zero-gravity boiling and condensation effects are handled by the com- 
bination of the THERM model and the liquid/ullage coupling model, for each of the 
THERM nodes. 

A basic block diagram of the zero-gravity computer model is shown in Fig. 3-9, 
with the interactions of THERM and the liquid/ullage coupling model solution sub- 
routine, ULLAGE 

An energy balance in the liquid gives: 

A AM + C sB AT s = AQ gu +AQ st + AQ gB 
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Fig. 3-9 Block Diagram for the Zero-Gravity Computer Models 
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where the three terms on the right represent respectively, the energy transfer at 
the liquid-vapor interface, at the wall, and in the bulk liquid, as determined by 
THERM. 

The bulk capacity, Cg^ = rCjJ «p 

,,t> i L *J L = T s 

From an energy balance on the total ullage volume V , we get: 

(C„T - C'T ) AM + C MAT + P AV = AQ 
v v p s' v sum 

It is assumed that the vapor behaves as a perfect gas over the limited pressures of 
interest. Thus the equation of state is used to relate tank pressure and vapor mass. 
The expanded form gives the following expression. 

ZRTAM + ZRMAT - V K c AT C - t> AV =0 

u S S S u 

dPg 

Kg in the above equation is the slope of the saturation line, i.e. , Kg = — . 

S 

The ullage volume is related to the liquid volume and the constant total volume, such 
that V = Vm - V T . When this is expanded in finite difference form, we have: 

U 1 Jj 


- v AM + AV = V T AM_ t -AM t AV 
L u L CL L L 

The four equations given above form a B matrix. The computer solves the matrix 
for the changes in total vapor mass, M , ullage temperature, T , saturation tem- 
perature, T , and pressure, P . 
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3. 3 FREE CONVECTION MODEL DESCRIPTION 


Figure 3-10 is a block diagram illustrating the interaction between the four subroutines 
within the free convection program. The problem is divided into a heating and cool- 
down phase. Transition between the two is set by the maximum and minimum pressures 
specified for tank pressure control. 

The basic approach is to determine the interface configuration and associated liquid 
height at the wall, for a given g level, tank size and ullage volume. This is then 
used to determine a psuedo- geometry for a flat bottomed, cylindrical tank that will 
yield the same modified Grashof (Gr*) number, liquid volume, and wall heat rates as the 
actual tank. This geometry conversion is accomplished in a GEOM subrountine. This 
geometry is then used in the STRAT subroutine to compute the pressure rise during 
the heating portion and in the COOLIT subroutine during the cooling phase. The 
ULLAGE subroutine is the same as that in the conduction model. 

3.3.1 GEOM Subroutine 

For any tank with hemispherical ends the volume is given by: 

V T = V/ U ° S [ 3 < L/ D> - X ] 

The ullage and liquid volumes are: 
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Fig. 3-10 Block Diagram for the Free Convection Computer Program 
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The volume and height relationship of a spherical segment is, 


h 


R + 2R cos 




, for k = 0, 1, 2 


where 


R = D/2 
= cos 




4 7T 



The correct height h of the spherical segment is the one of the three values in the 
range 0 < h^< R . 

The liquid elevation Z of the flat interface is found in the following manner for any 
tank fill level. 


If 


If 


If 


V L < 2/3 ttR 3 , then Z = h 

2/3 ttR 3 < V L < 2/3 ttR 3 [3 (L/D) - 2^ then Z 

V L >2/3 ir R 3 [3 (L/D) - 2 ] then set V g = V u 
Z = L - h 


= R + 


2/3 R ^3 (L/D) - 2 J 


in the spherical segment solution 


The interface shape due to low bond number accelerations must also be defined. This 
is accomplished with the ellipsoidal shape approximation of Satterlee (Ref. 9) from 
which the rise, at the wall, above the flat interface is found. For the propellants LH^ 
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and L0 9 , the contact angle 6« 0°, and this yields a very simple ellipsoidal 
approximation of the interface shape. 



o 


For 0=0° in a cylindrical section, then a = 1.0 and cos 0 — 1.0, and Ref. 9 
gives the following expressions: 



£ 3 + £ - (B <7 + 2) = 0, which has the form 


+ a £ + b = 0, where a = 1 and b — - (B<j + 2) 
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( BO + 2 ) 
2 


♦ v®^p> - 


(Bo + 2 ) . X 

4 27 


1 

p = — , where p = semi-minor axis of the ellipse, 


Additional length AZ along the wall above the flat interface is; 


= 2/3(3, or AZ - 2/3 p R 

XV •» 

In cylindrical section or the bottom hemisphere tjiis would also be increased elevation. 
However, in the upper hemisphere, the following relation for A Z is used where h< R : 

A Z = 2/3 p h 


This always insures some dry tank wall at the ullage, which would effectively be present 
with thin liquid layers. 

The liquid/vapor interface area is based upon the computed p value and the area of the 
flat liquid/vapor interface A^, . In the cylindrical section, the flat surface radius, 

Rp is: 


R f = D/2 


In the hemispherical ends, the radius of the flat interface R^, is 

R F = \Jr 2 - (R - h) 2 = y 2 Rh - h 2 

The flat interface area A_. is then: 

F 
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The actual interface area A is: 


A = (1+/3)A f 

The tank wall areas for bottom heating Ag, side wall heating A w and dry tank wall 
ullage A . can now be defined. If 6 g is the angle from the tank centerline at which 
the boundary layer starts, then this height hg above the tank bottom is: 
hg = R(1 - cos fl B ) 


The three areas of the 'tank wall are then: 


A B 7rDlL B 


A L a b + A w 


A = 7T D 
u 


A W = * D 


L - (Z+AZ) 


= * d [ d (f)- 


(Z+AZ) 


(Z+AZ) 




A t = ttDL = A b + A w + A u - ttD 


hg + L - (Z + AZ) + (Z+AZ) - hg 


= ttDL 


The boundary layer run length x along the wall can now be computed, depending on the 
value of (Z+AZ). 


If (Z+AZ) > R and (Z+AZ) < (L - R), then: 

X = R (ir/2 - 0g) + (Z+AZ) - R 

A different approach is required if the interface at the tank wall is located in either 
hemispherical end. 


If (Z+AZ) < R, then: 


Q = COS 
u 


-1 


(Z+AZ) 

R 


X - R <* u " *B> 
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If (Z+AZ) > (L - R) 

L - (Z+AZ) 

R 

X - R (tt/ 2 - 0 B ) + (L- D) + R0 u 

The above relations complete the geometry considerations required for the actual 
tank . 


$ = sm 

u 


-1 | (Z+AZ) - (L-R) 
R 


= sm 


-1 


3.3.2 STRAT Subroutine 

a 

The heating portion of the cycle, uses the Tellep and Harper Stratification Model I 
(Ref. 10). The actual tank is transformed to a cylindrical tank with a flat bottom, 
that will yield the same modified Grashof No. , Gr*, same liquid volume, , and 
same wall heat rates, q w , and q^ as the actual tank. With a uniform wall heat rate 
of on the actual tank, the parameters to be held constant from the actual tank are: 

* g*gp <C x4 P 

“ kv 2 


where; 


v = 


= liquid volume in the tank, which is kept current with time 

a s a M A 
TV vV w 

q B = ^ A B 

For the pure cylindrical tank (subscript c), these parameters can be expressed as: 


x 4 

C 


„ T- 

Gr - 

c 


kv 2 



V T = tt/ 4 D 2 
L c 

X 

c 


<* 

ii 

= q ^ * 

ttD X 
c c 

c c 

c 


q B = q B A B 

11 

bf- 

• 

tt/ 4 D 2 
c 
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The first three relationships above must be solved simultaneously for , D c and 

c 

X c for the pure cylindrical tank model, and yields: 


X 

c 


/kv 2 Gr *' W 3 - 5 
\ g+g o' iq w/ 



D 


c 



and 




o * 

Gr 

* 4 

* * 0 pX c 


The bottom heat flux q" can be then fixed as: 

C 



C 




The Model I stratification solutions of Tellep and Harper have both laminar and 
turbulent free convection mechanisms. However, it is assumed that translation 
occurs at one value of the Rayleigh number, where 


Ra 


* 


Gr *Pr 



for all liquid properties, then there would be a sharp discontinuity between the 

* 

laminar and turbulent solutions at the transition Ra number. To get around this 

* 

problem, transition is allowed to occur over a Rayleigh number range of Ra 
10 1:L < Ra < Ra MAX = 10 13 . For Ra values below this range, the pure laminar 
solution is used, whereas for values above this range, the pure turbulent solution 
is used. For solutions in this range, the parameter K is averaged in a weighted 
manner with respect to the Rayleigh number: 
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K 


ave 


( Ra MAX " Ra ) 

1 k lam + • 

/ * 

Ra - . 

Ra MIN ) 

1 k turb 

i 

( Ra MAX ~ Ra Mm) 

1 


In the Tellep and Harper Model I, the solution is divided into two regimes, depending 
of the mass startified layer A, has reached the bottom of the tank or not. The time 
parameter for these solutions is: 



The 


Ai 

X 


solution before 


x 

c 


1.0 has the following solution: 


| 

X c./LAM 


= 1 - 


1 - 0.630 



0 


Gr* ^ 

„ 0.388 
Pr 


5 



TURB 


1 - 


1 + 0.0924 



\ 

Gr 

1 + 0.443 Pr 2//3 


2/7 


-7 


The time <f > ^ required for the mass layer to reach bottom has been solved from the 
above equations: 



LAM 


Pr 


0.388 


1.5873 



Gr 


* 1/5 


L\ _5.75Pr 2/3 

1 + 0.443 Pr 2 ^ 3 

\ A) TURB /X\ 

* 

Gr J 

W 
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For 0 < 0 ^ 


Al 


X 


= 1 - 


c/LAM 


1 - 0.630 


X 

< 

R 


Gr V* 


Pr 


0.388 



With the appropriate value of 




, expressions for the three desired quantities 


T 0 , A 0 and follow. T is the initial liquid temperature at time t = 0 . 
S 2 13 o 
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When the stratified layer reaches the bottom of the tank (4> > 4^) an Y additional heat 
is asumed to be distributed uniformly, thus maintaining the shape of the temperature 
profile that occurs at time 4 = 4 ^* 

Below the temperature stratified layer A the liquid temperature is a uniform T^. 
In the stratified layer, the temperature distribution will be: 



where y* = —r ~ — , y is the depth below the liquid surface, and the equation applies 
A 2 

0< y*< 1 . 

Using the Tellep and Harper expression for f (y*) yields: 




e 

(1 - y*) where e = exponent. 


The exponent e is defined by fixing the value of the energy integral F , and using the 
relation: 


/I 

f(y*) dy* = f 

n 


1 e _n-v*i e+1 

(1 - y*) dy* = — * ^ 


e + 1 
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F 


+ J1L 

e + 1 


e. + 1 


L_ 

e + 1 


F ( e + 1) = 1 

Fe = 1 - F 


This completes the model of temperature stratification due to tank wall heating in the 
equivalent cylindrical tank with no cooling present. 

3.3.3 COO LIT Subroutine 

During cooldown, cold liquid stratification will occur to the bottom of the tank, and 
bulk boiling will occur in the upper warm stratified layer. The bulk boiling near the 
top is treated by the liquid ullage coupling model. For a constant cooling tube tempera- 
ture T and a decreasing mean liquid temperature, T t , the mean wall cooling rate 
will decrease with time. Since the cooling model must always result in net cooling 
to the liquid, a test must be made to insure this effect. Based upon a unit area of 
liquid/tank wall surface area, where h - p 1 - and h T „. = f (Ra) for the free 

CW -tv -LiW 

CW 

convection, then the requirement for net liquid cooling would be: 


h 


cw 



L min 



Where T^ = saturation temperature at minimum pressure. 

An energy balance on a unit area of liquid tank wall surface area would be as follows, 
assuming T^ = mean liquid temp in the tank. 

h <T„_ - T ) = q" + h T „,(T T - T n ,), and 
cw W c Tv LW L W 
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(T - T ) = 
' L w' 


h (T_ - T ) - q M 
cw x L c ; Av 

(h + h TTT .) 

' cw LW' 

The free convection heat transfer coefficient h LW is based upon correlations for 
constant AT. The actual tank dimensions are used to compute Ra and the tank 
cooling: 

s* s 0 X 3 

Gr g ( t l " T w ) = regular Grashof number based upon AT. 

/“i i j 

W ** 

Pr = — — = liquid Prandtl number. 

Ra = Gr Pr = regular Rayleigh number based on AT. 

Based upon the Rayleigh number magnitude, the Nusselt number is calculated from the 
following: 

Nu = — ^ — - 0. 105 Ra^^ For Ra > 10 ^ (Turbulent) 

1/4 4 q 

Nu = 0.59 Ra For 10 < R< 10 (Laminar) 

4 

Nu = 5. 90 For Ra< 10 (Viscous and Conduction) 

The value of h^^y is updated for each computed time step based upon the initial 
value of (T^ - T w ). Since the wall temperature is not known at the start of the 
cooling cycle, the following method is used to estimate it: 

1* Set h^y = 0 and estimate first value of (T - T ). 

2. Use this (T L - T w ) to compute Ra, Nu, and h LW . 

3. Use this h^^y to compute next the value of (T^ - T ). 

After this first cooling iteration, the current value of (T_ - T ) is used to predict a 

L W 

value of for the next time step. 
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The energy transferred from the total liquid mass is. 


q L " h LW A L < T L _ T w) 


This mean drop in temperature of the liquid T is found from q^ 

dT dT q 

c l TE” = q L> hence TT- = ~cl 


AT. 


dT } 

dt 


, ^At,j whereAt - compute time step. 


This AT l is successively added to the liquid temperature distribution. During cooling, 
a deeper layer of liquid at the top will be at saturation temperature T g , which will 
have the capacity C gB . For the matrix solution in the liquid-ullage coupling sub- 
routine, the average heat rate q gB to this layer must be fixed from the overall cooling 
value q^. 



The mean liquid temperature T^ is found during heating and cooling in the following 
manner. The temperatures of the upper liquid layers, the remaining stratified layer, 
and the bulk liquid at the tank bottom are averaged, with respect to capacity, to yield 
temperature T^. A running average is kept for this quantity. During cooling the 
liquid withdrawn for the cooling tubes is assumed to come from the bottom of the tank 
at this mean bulk liquid temperature. 


This program computes the liquid and tank geometry only once, at the minium pressure, 
during each cycle. As with the conduction model, the liquid capacity is also updated 
only at each restart of a pressure cycle. 
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3.4 ANALYTICAL PRESSURE HISTORIES 


The previously described conduction and convection models have been used to 
generate pressure and temperature histories that span the range of tank sizes, orbital 
heat fluxes and gravity levels of interest. These are presented to illustrate the in- 
fluence of the various parameters. 

3.4.1 Conduction Results 

Results are presented for the tanks with the largest and smallest capacity, and for 
both 5 percent and 90 percent ullages. Tanks sizes are, respectively: D = 6.86 M 
(22.5 ft) xL/D - 4, and D = 1. 22 M (4 ft) x L/D = 1. The nodal structures are 
shown in Figs. 3-11 through 3-14. 

Figures 3-15, a through f, show the resulting pressure and temperature histories in 

hydrogen for the large tank configuration shown on Fig. 3-11. Figures a through c are 

2 2 

for a heating rate of 3. 15 w/m" (1.0 Btu/hr-ft ), and cooling resistances of .09, .35, 

. 88 m2 ° - K (0.5, 2.0, and 5.0 . 

w ' Btu 

2 

Figures d through f are for similar resistances but the heating rate is 0. 1 Btu/hr-ft . 
For all the runs, the operating pressure was (15-18 psi). 

It can be seen from Fig. 3-15 a through c that the heat exchanger (R^,^) does affect the 
vent cycle period, and in fact, for the largest resistance (wide tube spacing) the cooling 
rates were insufficient to pull the pressure back to the starting pressure (i.e. , could 
not control). 

The liquid node temperatures are referenced to Fig. 3-11. They contain temperatures 
along a radius from the wall to the tank centerline and axially from the bottom to the 
liquid-ullage interface. They confirm that the changes are limited to the region near the 
tank wall. They also indicate that subcooling is more severe for the first cycle. This 
appears to be a consequence of the starting condition wherein all vapor is contained in 
the ullage. However, part of the initial ullage condenses and vapor is generated along 

3-34 


LOCKHEED MISSILES & SPACE COMPANY 




















Fig. 3-12 Liquid Node Structure For Conduction Model With a Cylindrical 
Tank - L/D = 4; 90% Ullage 
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Fig. 3-13 liquid Node Structure For Conduction Model With a Spherical 
Tank; 5% Ullage 
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ULLAGE NODE 



Fig. 3-14 Liquid Node Structure For Conduction Model With a Spherical 
Tank: 90% Ullage 
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the wall during the first cycle. Some subsurface vapor exists throughout the ensuing 
cycles. The subcooling is reflected on the pressure history curves as a lag in the 
start of a subsequent rise in pressure following cessation of cooling which occurs as 
soon as the pressure drops to minimum level. During this lag time the heat coming 
into the tank is going into warming up the subcooled field. 

The results shown in Figs. 3-16 a and f are for the same tank with 90 percent ullage. 

With 5 percent ullage, the vapor constitutes approximately 0. 1 percent of the total 

thermal capacity whereas it represents 17 percent for the large ullage. The time 

period for a vent cycle is indicative of the extent to which energy is being distributed 

in the liquid. Uniform distribution would give much longer periods for the small 

ullage (large liquid capacitance). However, when the energy exchange is confined to 

a small portion of liquid, phase change and compression (or expansion) of the ullage 

results in rates of pressure change which are more nearly proportional to the ullage 

volume. In this situation larger vent cycle time occurs for the high ullage volume. 

The latter situation is found to be the case when the heat flux is 3. 15 w/m (1 Btu/ 

2 

Hr-Ft ), for which the vent periods were approximately an order of magnitude greater 

with 90 percent ullage than with 5 percent ullage. At the lower heat flux (0. 1 Btu/ 

2 

Hr- Ft ), the vent cycle is only 20 percent longer for the large ullage. The implication 
of this is that non-uniform (stratification and subcooling is less severe at the high 
ullage volumes. This is substantiated by the temperature histories. 

The temperature histories are shown in Figs, 3-17 and 3-18 for the small tank. The 
trends are what one might expect relative to the large tank results. The thermal 
oscillations penetrate further into the tank. In particular, notice that at the lower heat 
rate, (Figs. 3-17 d through f) even the center most liquid temperature (node 39) is 
affected. It continues to warm up for a few hours after the cooling is started until 
such time as the tank pressure drops to the point where boiling takes place (on Fig. I7f 
this occurs at 59 hours). During the cooldown, subcooling occurs in the three nodes 
nearest the wall (nodes 14, 23, and 27). There is also excess cooling of the ullage 
vapor which produced a pressure undershoot to 10 n/m (14.6 psi). However, as with 
the large tank, the excess cooling was recouped by absorbing additional heat in the 
prolonged warmup portion of the next cycle. 

3-39 


LOCKHEED MISSILES & SPACE COMPANY 



For the purpose of comparison, a few selected cases are shown for oxy.gen on Figs. 3-19 
through 3-22. For the small tank all liquid nodes have essentially the same temperature 
history, with 90 percent ullage, indicating a virtual absence of stratification. 


The liquid oxygen runs display characteristics similar to those found with hydrogen. 

The most noticeable difference is that a higher resistance (wider spacing) heat exchanger 
design can be effectively used with oxygen. With hydrogen, a design resistance of 

rn 20xf Q rt 

0. 88 (5 hr-ft - R/Btu) was so high that the heat removal rate was insufficient 


+ a n-PfciA+i xml tt nnntvAl fKci foMlr nr'occni'o 

V \J ViiVCVl Wlj WJL Vliv VUililV VMWIAJU v ( 


With rw\Tcrc*r\ tntiil tAmnArntnrp 


between the bulk fluid temperature and that of the expanded fluid in the coolant tubes 
is 9R° compared to 4.4R 0 for hydrogen. This, coupled with the lower thermal diffusivity 
of oxygen indicates that the heat exchanger wall resistance can take on a higher value 
for oxygen than for hydrogen before it becomes the limiting factor in pressure control. 


The hydrogen data indicated that an order of magnitude decrease in heat flux (from 

9 

3 to 0. 3 w/m“) increased the cycle time nearly two orders of magnitude for the 5 per- 
cent ullage. A similar effect was obtained with oxygen. The cycle time for liquid 
oxygen is approximately 6 to 10 times that for hydrogen with the full tank. The 
thermal capacity of LOg, for a 3 psi control band, is 8. 5 times that for hydrogen. 
With the high ullage, the cycle time is increased by a factor of 10 for a similar de- 
crease in heat flux. This is consistent with the hydrogen runs. 


The venting inefficiency at the end of a cooling period is represented by the amount of 
excess mass that was removed from the tank to produce the subcooled condition. This 
excess vented mass can be expressed as a percentage of that which would be 
vented up to that point in time for a completely mixed propellant? i.e. , 
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where 


inefficiency factor 
Excess vented propellant 
Vented propellant for mixed system 
Amount of subcooling 

Total heat input that would produce w yQ ' rv 

Figures 3-23 and 3-24 show this inefficiency factor as a function of time. It decays 
exponentially, as one might expect, primarily due to linear increase in the denominator 
(Qrp). It is interesting to note that the effect of tends to converge as time pro- 

gresses. This too, is because of the fact that subcooling becomes less important as the 
mission proceeds, and consequently, perturbations due to heat exchanger design be- 
come second order effects. 

Figure 3-25 presents the same data correlated against vent cycle number. This tends 
to pull the data for the two heat rates together in spite of the fact that vent period and 
frequency differ by two orders of magnitude. This indicates that the performance is 
relatively independent of frequency, as long as the design is adequate to provide control 
(e.g. , if q = 1. 0 Btu/hr-ft 2 , R cw 5). Moreover, a space mission requiring this 
size tank, and these heat rates would likely be of long duration and require considerably 
more than 5 vent cycles which implies that a passive heat exchanger will not result in 
a significant weight penalty, even in the absence of convection. Also, it is important 
to note that the subcooling penalty is computed at the end of a cooling period and does, 
in fact, represent a penalty only if the mission were completed at that time. If the 
mission is not completed, the non-vent portion of the next cycle is prolonged by virtue 
of the increased thermal capacitance of the subcooled fluid. Thus, the weight penalty 
is a discontinuous function which is negated in the pressure rise portion of each pres- 
sure cycle. Inasmuch as these conduction-dominated cases represent the most severe 
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conditions for operation of a wall-mounted heat exchanger, there is no doubt that this 
concept can be effectively used in a thermodynamic vent system for hydrogen and 
oxygen propellant tanks that are to spend prolonged period in space. It is also applic- 
able to non-cryogenic fluids, but the high vapor pressure cryogens present the most 
difficulty in storage and venting during low gravity space flight. 


3.4. 2 Convection Results 


Figures 3-26 through 3-30 illustrate the type of results obtained with the convection 
model. These are for a tank that is 22. 5 feet in diamter and 90 feet tall. Figures 3-26 
through 3-28 indicate the relative effect of the thermal resistance in the tank wall 
(R ) between cooling tube attachment points. 


A large resistance indicates wider spacing between tubes, fewer tubes, and thus lower 
heat exchanger weight. As expected, a larger resistance, required a longer time to 
pull the pressure down. Also, it can be seen that temperatures within the body of 
the liquid tend to decrease throughout the mission, and this subcooling is reduced with 
increasing values of R cw * 


Figures 3-29 and 3-30, in conjunction with Fig, 3-27 illustrate that the increased 

stratification associated with higher g-levels tends to shorten the pressure cycle as 

well as the amount of liquid subcooling. The average cooling flux is shown on Fig. 3-31, 

as a function of g-Ievel for a design having a resistance value of R = 0. 17 m ^°^- 
ft2-OR cw w 


(2 


Btu 


). 


Also shown is the cooling flux determined with the zero-g conduction model for the 

same tank and heat exchanger. The indication is that the latter model is appropriate 
“6 

below 10 standard gravities. 
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Fig. 3-15 Pressure and Temperature Characteristics in 
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3.5 HEAT EXCHANGER DESIGN MODEL 


The conduction and convection models are basically performance analysis tools that 

describe the processes in the tank fluids for a given heat exchanger parameter (R cw ). 

R is the thermal resistance between the tank wall and the vent fluid in the heat 
cw 

exchanger as illustrated on the sketch below. 

i/hi , 

9 -WSA # vvSSv • 

T L T w T c 

The heat exchanger geometry (tube spacing and diameter) must be derived from con- 
siderations of conduction in the tank wall and the pressure drop in the coolant tubes. 

To be effective, any thermodynamic vent system must be designed to avoid excessive 
local heat transfer which might vaporize all of the refrigerant and inhibit effective cool- 
ing and depressurization of the ullage. With an active system, such as that developed 
in contracts NAS 3-7942 and NAS 3-12033, this is accomplished with use of a mixer in 
the propellant. For the tank wall exchanger the tubes will be attached to the tank wall 
at discreet points to provide uniform heat extraction over the tank surface. This divides 
the tank wall into a large number of fins. Heat is convected or conducted into the tank 
wall and then conducted radially to the attachment point. The heat to each attachment 
is given by: 


K,(a V/3) 

Q p - 2*kt (a < T c- T L> 


The total heat rate is: 


Q t = NQp = (hA g rj> (T c - T l ) 
These equations are used to define an effectiveness which is: 


V 


_N 

A 

s 


2 7rkt (a/ /3) K Q (a / f3) 
h 
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It can be shown that the Bessel function component of the above expression for the 
range of g levels and tank thicknesses of interest in this study is relatively constant 
at 0.5; i. e. , 


K ± (a//3) 

K (a/B) 
0 b 


(a/0) = 0.5 for a//3 <0.7 


If the attachment spacing, S , is small compared to the basic dimension of the tank 

2 

wall, then one can approximate the tank surface area as NS . Making this approxima- 


tion the effectiveness for the externally mounted heat exchanger thus reduces to: 


V “ (a/S) X 77 


For an internal heat exchanger heat will be transferred to the wall and also directly to 
the tubes. If it is assumed that there is no synergistic effects, and assuming that the 
tube dimension and tank attachment radius are the same, the effectiveness for an 
internal exchanger becomes: 

t) 2 = ^ + 2tt (a/S) 


Figure 3-32 gives a graphical comparison of these two mounting arrangements. It indi- 
cates that on a theoretical basis, the internal mounting results in fewer tube attachments 
for a given effectiveness. Operationally, however, this concept is susceptible to fail- 
ures if a leak develops. On Contract NAS 3-12033 leaks developed between the expansion 
unit and the heat exchanger which resulted in by passing the former and eliminating the 
temperature difference needed for the exchanger to operate. The tank pressure could 
not be controlled when that happened. In contrast, a leak in the external system will 
result in local degradation of the insulation system but will allow a mission to proceed 
at a reduced capability. 

The difference between the internal and external effectiveness indicates that most of the 
heat is transferred directly to the tubes rather than to the tank wall. This tends to 
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negate the purpose of the discreet point attachment concept for assuring uniform 
energy extraction. 


Figure 3-32 provides a relationship between wall fin effectiveness and the geometry 

ratio (a/S) . The effectiveness is also related to the wall resistance R by: 

cw 


rj = 1 


R q T 

cw L 

(T t - T 
' L c 


) 


The average cooling flux, q T , and temperatures were obtained from the computer 
results with the zero-gravity model, and the effectiveness computed with the above 
equation. Figure 3-33 shows the effectiveness obtained in this manner, using results 
for the 4 foot sphere and the large tank (22.5 x 90 feet). The indicated band repre- 
sents scatter in the data from one pressure cycle to the other through the course of 
mission history. Figure 3-34 is a cross plot from Figures 3-32 and 3-33 to relate R 
to the tube radius to spacing ratio. 


It might be noted that the intercept (O^R ,) on Figure 3-33 corresponds to a heat 

flux conducted to the coolant that just equals the imposed heat flux from external 

environment. Therefore, any design resistance greater than the intercept value 

would result in a net heat into the propellant, even during the cooling portion of a 

pressure cycle. There would be a concomitant inability to control the tank pressure. 

2 

This, indeed, was the case when the imposed heat flux was 1 Btu/hr-ft and the 

o 2 

value of R was 5 R-ft -hr/Btu for the hydrogen tanks. 

R cw is increased by increasing the spacing between attachments. Consequently, the 
number of tube coils and tube weight decrease with increasing R cw * which indicates 
the desirability of operating near the intercept value and the attendant continuous 
venting mode. 

In reality, continuous operation is impractical with a passive system because it does 
not provide any flexibility to accommodate changes in the environmental heating. For 
example, if a system is designed to maintain continuous, constant pressure with a 
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ank Wall Resistance, R cw , -°K-m 2 / 



Tank Wall Resistance,R cw> -°K-m 



0 CM 0.2 0.3 O.li 0.5 
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Figure 3-34 Effect of, Attachment Radius to Spacing Ratio on Tank Wall 
Resistance — External Tubes 
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nominal heat flux, it will lose control if the actual heating during any portion of the 
mission is higher than the nominal value. Consequently, one will design for con- 
tinuous operation at the highest potential heat flux and except intermittant operation 
at nominal anticipated values. 

Another consideration in selecting a design is the practical limitation of other com- 
ponents in the thermodynamic vent system (i. e. , expansion unit and flow control 
values). During the conduct of NAS 3-7942 it was estimated that design flow rates 
"below 1 pound per hour were not realistic from the standpoint of automatic operation. 

To flow continuously at this rate corresponds to an effective heat rate of 190 Btu/hr. 
Thus, for an actual heating flux of 3.15 w/m 2 (1 Btu/hr-ft 2 ), continuous venting would 
be practical only for tanks having greater than 190 sq. feet of surface area. Conversely, 
continuous venting (effectively) of smaller tanks would be realistic only if the actual 
heating flux were greater than 3. 15 w/m 2 . Figure 3-35 incorporates these consider- 
ations into a relationship between maximum R and minimum tank surface area for 

cw 

which continuous operation might be feasible. For a given tank, the value of R max 

cw 

obtained from Fig. 3-35 can be compared to the intercept value on Fig. 3-33 at the 
appropriate heat flux. 


The smaller of the two values of R cw is the appropriate design value. For example, 
the surface area of the 1.2 meter (4 ft) sphere is 4. 7 m 2 (50 sq. feet). From Fig. 3-35 
R^ max is 0. 17. If the maximum anticipated heat flux is 3. 15, then Fig. 3-33 gives 

an intercept value of 0. 6 . Thus, R = 0. 17 is appropriate design value and the 

cw 9 9 

system would operate intermittently for a heat flux of 3. 15 w/m (1 Btu/hr-ft ) . 

Figure 3-34 provides the correlation between R cw and a/S . The ratio between tube 
weight and tank membrane weight can be expressed as: 



x 2 tt aStj. 


q 2 

p x S 
H x 


t 

s 



(a/S) 
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Figure 3-35 Relationship Between Minimum Tank Surface Area and Maximum R cw For 
Continuous Operation 








This weight ratio is shown in Fig. 3-36, for aluminum tubing on aluminum tanks. 
The magnitude of the heat exchanger tubing weight is given on Fig. 3-37 for specific 
tanks. The tube wall thickness used for Fig. 3-37 was 6 mils. 

We must look to the pressure drop to actually size the tube. The pressure drop is 
comprised of three elements: That due to friction; due to gravity; and that due to 
acceleration when the fluid is being vaporized. The expression which describes the 
pressure drop from inlet to exit of the tube is derived in Appendix B. It is: 


[\ 

2 

[ f i G v i 

1 + -lj 

, G 2 ve/ V e \ dx _ *V tae ] 

4 A S c 

g c \v 1 / d l T Ve 1 +X (VgA^-l) 




_2 dv 

I 



1 + X — — - = t 

f 

0 

- 


g c dP J 


The solution can only be obtained by repeated iterations for specific mass fluxes and 
assumed geometries. However, the bracketed quantity in the denominator is the 
square of the Mach number in the vapor and must be less than 1 for a constant diam- 
eter tube to avoid choking. This fact can be used to estimate the absolute minimum 
tube size to be considered for a given tank. This would constitute a logical starting 
point for iteration on the complete solution in a very specific application. 


Since the mass flux is the vent rate divided by the tube flow area, and the vent rate 
is determined by the thermal design parameters, one can derive the following ex- 
pression for the tube radius. 


V7T 

s 


< T L “ T c> (1 " 


1/2 


G 


7T 


R 


cw 


When G is large enough to choke at the exit, the tube radius is the minimum possible. 
The larger the pressure expansion between the tank and the heat exchanger inlet, the 
larger will be the temperature differential, but the smaller will be the possible pres- 
sure drop and mass flux. Assuming the expansion to 4 psi, which is the same as that 
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used for the active system in contract NAS 3-7942, one obtains the relationship be- 
tween the normalized tube radius (normalized with tank surface area) and the tank 
wall design resistance that is shown on Fig. 3-38. For a given tank, the minimum 
tube spacing and radius can then be determined as shown on Figure 3-39. 
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Figure 3-37 Heat Exchanger Tube Weight Variation With R , Material: 
Aluminum, 6 mils cw 


3-125 


LOCKHEED MISSILES & SPACE COMPANY. INC. 



LOCKHEED MISSILES & SPACE COMPANY, INC. 


co 

i 

h- 1 

N) 

Cl 



0.005 


0 . 001 + 


0.003 


0.002 


0.001 


0 



R cw - °R-ft?-hr/Btu 



Figure 3-38 Normalized Heat Exchanger Tube Radius as a Function of Tank Wall 
Resistance 
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Fig. 3-39 Minimum Tube Radius and Spacing as a Function of Wall 
Resistance and Tank Size 
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Section 4 

EXPERIMENTAL PROGRAM 


This test program was designed to provide data needed to verify or modify the analytical 
models. Basically, this required determination of cooling rates, temperature dis- 
tributions, and pressure responses with various input conditions selected to simulate 
a range of low gravity environments. 

Testing was performed with three non-cryogenic and two cryogenic fluids to cover 
approximately a four decade range in Rayleigh number. Both a cylindrical and a 
spherical test apparatus were used and tests were conducted at various ullage volumes 
to determine the extent to which performance is effected by geometry. This section 
of the report contains a description of the test apparatus; instrumentation, conditions, 
procedures and results. 

4.1 NON-CRYOGENIC TESTS 

The three fluids used for these tests were: Freon-11; water; and 30 weight motor oil. 
The fluid properties give approximately 3 orders of magnitude variation in Rayleigh 
number. An attempt was made to extend the range by packing polyethelene spheres in 
the tanks. Several papers (e. g. , Ref. 14) have shown that the onset of convection can 
be retarded by introducing a porous media, for which the effective Rayleigh number is 

t 

defined as: 


Ra = Ra 
eff o 


..Li i «i/L )' 

Ll50 (l-e) 2 


"1 

■J 


where: 


Ra Q is the basic fluid Rayleigh 
6 is media porosity (0.476 for packed spheres) 
d is diameter of packing spheres 
L is height of sphere packing 
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Both 3/4 and 1-1/2 inch diameter polyethelene spheres were used. Table 4-1 indicates 
the range of Rayleigh Nos. that were anticipated with these fluids, two tanks, and two 
porous media sizes. 


Table 4-1 

ESTIMATED RAYLEIGH NOS. FOR COOLING 


Tank — 


Cylinder 


Sphere 

Fluid — 

Water 

Oil 

Freon 

Freon 

Media , 
Size 

None 

2xl0 10 

2xl0 8 

5xl0 U 

GxlO 11 

3/4 

4 

5x10 

4xl0 2 

2xl0 6 

2xl0 6 

1-1/2 

2xl0 5 

2xl0 3 

6xl0 6 

7xl0 6 


Table 4-2 presents the matrix of test configurations that were included in these non- 
cryogenic tests. 


Table 4-2 


NON-CRYOGENIC TEST MATRIX 


Tank — 


Cylinder 


Sphere 

Fluid -*■ 

Water 

Oil 

Freon 

Freon 

Media , 
Size * 

None 

3/4 

1-1/2 

1,2,3 . 
1,2,3 

1,2,3 

1,2,3 

1,2, 3, 4 
1,2 

1,2, 3,4 

1,2,4 

1,2,4 


(a) Numbers in the matrix columns refer to the follow- 
ing ullage configurations 

1 '5% ullage, Flat Interface 

2 10% ullage, Flat Interface 

3 10% ullage, Spherical Interface 

4 50% ullage, Flat Interface 
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4.1.1 Non-Cryogenie Apparatus 


The two test articles employed in these tests are shown in Fig. 4-1 and 4-2. The 
first is a right circular cylinder made from 12 inch stainless pipe. The wall thickness 
is 0.156 inches. The heat exchanger is made from l/2 inch copper tubing. It con- 
sists of 7 coils attached to the tank with clips equally spaced at 3.5 inches. A niehrome 
wire heater is uniformly positioned between the heat exchanger coils to rapidly heat 
up the tank contents. 

As illustrated in Fig. 4-1, the cylindrical, apparatus had an 8 inch diameter aluminum 
sphere installed for part of the tests to simulate a submerged ullage bubble which was 
cooled indirectly by the heat exchanger. This sphere constituted 10 percent of the tank 
volume. The bubble had a neck made from 1/4 inch tubing which penetrated the flange 
cover in the top of the tank. This neck provided the physical support for the sphere 
and the fluid connection between it and the main tank. An isolation valve (IV01) was 
placed in the connecting line between the two tanks so that they could be filled with 
a common fill line and drained individually. For a spherical ullage simulation the 
tank was filled to the top. Drain valve DV01 was then opened, and when the sphere 
was empty, the isolation valve IV01 and DV01 were then closed to lockup the gas 
pocket. For flat ullage tests, the isolation valve was left open and the prescribed 
quantity of liquid was drained from the main tank through DV02. 

The major element of the spherical test apparatus is a 22.5 inch diameter hydrogen 
tank that had been previously developed for an ID program. It was made from 2219 
aluminum with a nominal membrane thickness of 0.019 inches. It had a 4 inch diam- 
eter opening with a mating, conoseal access cover. This prohibited the use of the 
8 inch spherical ullage in this tank. 

The heat exchanger consisted of 9 coils uniformly spaced, and attached to the tank at 
3.4 inch intervals. This provided approximately the same thermal resistance as for 
the cylindrical apparatus. 

Figure 4-3 is a schematic illustration of the test set-up for these non-cryogenic tests. 
The entire test apparatus was placed in an insulated box to minimize heating from the 
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Figure 4-1 Non-Cryogenic Cylindrical Test Apparatus 
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(t) Temperature measurement 
(g) Pressure measurement 

(a) Schematic 


(b) Hardware 


Figure U-2 Spherical Test Apparatus 
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(T) SIGNIFIES TEIPEiiATUjRE IlEASUREfENT 

Fig. 4-3 Non- Cryogenic Test Diagram 






external environment. The insulation also supported the test tanks. As illustrated, 
the heat sink (simulating vent fluid in the heat exchanger) was provided by water 
flowing in a recirculation loop through an air conditioning unit. The controller on 
this air conditioner permitted coolant temperature cycling between 40 and 48 degrees. 
This temperature was lowered and the variation was reduced by routing the coolant 
water through approximately 35 feet of copper tubing that was immersed in an ice bath. 


The location of the temperature and pressure measurements are shown to scale on 
Figs. 4-1 and 4-2. The tank wall was instrumented to give a temperature profile be- 
tween tube attachment points at two locations. A single temperature was also obtained 
between the other tube coils to give a measure of the variation from bottom to top of 
the tank. Liquid temperatures were measured at several points to give a good defini- 
tion of the profile and for establishment of the average temperature differential history 
during the cooldown. Table 4-3 is a listing of the instrumentation for these non- 
cryogenic tests. Figure 4-4 is a photograph of the internal instrument probe. 


Table 4-3 

NON-CRYOGENIC INSTRUMENTATION 


Item 

Type 

Range 

No. of 
Measures 

Elevation 

(inches) 

No. of 
Measures 

Elevation 

Cylinder 

Sphere 

Ullage 

Pressure 

C.E.C. 

0-25 psia 

1 

23 

1 

22 

Internal 

Fluid 

Temp, 

PRT Ch-Cn 
Differential 
Thermocouple 

30-540°F 

0-30°F 

1 . 

7 

| 

1 

7 

1 

5,10,9,11 

18,19,20 

Wall Temp 

Cu-Cu 

Thermocouples 

0-30°F 

10 

HI 

8 

. 5, 1, 1. 3, 

3.2,6,9,16 

19,20,21 

Coolant 
Temp In 

PRT 

30-540°F 

1 

— 

1 

— 

Coolant 

Temp 

Out 

PRT 

30-540°F 

1 


1 

< 

Coolant 

Flow 

Turbine Flow 
Meter 

0-4 gpm 

1 

— 

1 

- 

""" 
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Fig. 4-4 Internal Instrumentation Probe for the 0.56 m (22 inch) Diameter Test 
Applications 
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4.1.2 Non-Cryogenic Test/Procedures 


Prior to conducting any tests, functional checks were made on all instrumentation, 
heaters, etc. The test article was slowly filled with liquid. A prescribed quantity 
of liquid was then drained from the tank to set up the desired ullage configuration. 

A typical sequence of operations is as follows: 

(1) Apply a uniform heat rate until the temperature in the top of the tank 
increases 10°F» 

(2) Record temperatures and heat rates until the desired temperature rise is 
achieved. 

(3) Initiate coolant flow at approximately l/2 gpm and record all temperatures 
as function of time until surface temperature returns to original value. 

(4) Drain the tank to the next level and repeat steps (1) through (3). 

4.2 CRYOGENIC TESTS 

The cryogenic tests were conducted in the 22 inch diameter sphere, using both liquid 
nitrogen and hydrogen. These fluids extend the Rayleigh No. range another order of 
magnitude. However, because of their very high vapor pressure they allow a meas- 
urable pressure response which provides the best indication of mass and energy trans 
ger at the interface, as well as providing a check on the validity of the pressure 
control model. 

For the non-cryogenic tests, the coolant in the heat exchanger was supplied from an 
external, controlled temperature water supply. The cryogenic tests, were intended 
to demonstrate the automatic control characteristics of a tank wall heat exchanger 
system. Tank fluid was extracted, and expanded across a Joule -Thompson device 
to give a temperature drop. This subcooled fluid was then vented through the heat 
exchanger to cool the tank contents. The approximate Rayleigh Numbers that were 
anticipated for the cryogenic tests are shown in Table 4-4. It was intended to include 
the 1-1/2 inch polyethelene sphere media. However, the non-cryogenic tests proved 
conclusively that this media did not significantly reduce the strength of the convention 
boundary layer. Thus, the porous media tests were not conducted with cryogens. 
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Table 4-4 

ANTICIPATED COOLING RAYLEIGH NUMBERS - CRYOGENIC 


CRYOGEN 

ln 2 

LH 2 

Ullage — % 

15 

50 

15 

50 

Rayleigh No. 

io 12 

3 x 10 11 

2 x 10 12 

11 

7 x 10 


4.2.1 Cryogenic Test Apparatus 

The cryogenic test apparatus is the same as the spherical article described above, 
with one minor change. Whereas for the non-cryogenic tests the coolant was supplied 
by the external loop, for these tests it was extracted directly from the tank. There- 
fore, the coolant inlet was connected to the Joule Thompson valve as illustrated in 
Fig. 4-5. 

The apparatus was installed in the vacuum chamber as shown in Fig. 4-6. In as much 
as the cryogen had to be expanded to less than l/2 an atmosphere to get sufficient tem- 
perature differential, a vacuum pump was used downstream of the heat exchanger to 
suck the vent fluid through. Vent fluid flow was measured with a Hastings cryogenic 
thermal mass meter. The quality of the vent fluid was determined indirectly by use of 
a calibrated heater in the vent line downstream of the heat exchanger. This heater 
was used to bring the vent fluid to a superheated condition. The exit quality was then 
calculated from measured power, flow rate, and temperature rise across the heater. 

The instrumentation for the cryogenic tests will be the same as previously described 
for the sphere in the Task H tests, except for the flow measurement device. The 
Hastings thermal mass meter, which has a range of 0 to 2 pounds per hour of 
hydrogen, was used to measure the flow rate through the tank wall heat exchanger. 

A Rosemount thermal mass meter located in the tank vent line was used for boiloff 
determination. 
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4.2.2 Cryogenic Test Procedures 


The cryogenic tests were conducted using the pressure switch that was used with the 
active LIQUID PROPELLANT THERMAL CONDITIONING SYSTEM (Contract NAS 
3-12033). A typical sequence of operations is as follows: 

(1) Install test article as per Fig. 4-5. 

(2) Evacuate chamber and test apparatus. 

(3) Open facility shutoff valve and initiate fill until liquid level is at 5% ullage. 

(4) Allow tank pressure to stabilize and determine equilibrium heat rate. 

(5) Place system on automatic operation with the heat exchanger system. 

(6) Record all temperatures and pressures and vent flow rate at prescribed 
intervals during the cooldown portion of a pressure cycle. 

(7) Record temperatures and pressure at prescribed intervals during the 
pressure rise portion on the cycle. 

(8) Drain the liquid level to the next ullage point and repeat steps (4) through 

( 7 ). 

(9) When all three ullage volumes have been completed with LN 2 , drain and 
evacuate the test article. Then repeat steps (3) through (8) with LH 9 . 
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Fig. 4-6 Facility Test Schematic — Cryogenic Tests 




4.3 TEST RESULTS AND ANALYSIS 


4.3.1 Discussion of Results 

A summary of all the test conditions is presented in Table 4-5. The initial liquid 
temperature shown is the average of the measurements in the liquid at the start of 
cooling. The coolant temperature is the average between the heat exchanger inlet and 
outlet values. The heat flux was determined from the heater power measurements for 
the cylindrical tank and from boiloff measurements for the cryogenic tests. 

The temperature histories presented in Figs. 4-7 through 4-23 illustrate the typical 
performance of the tank wall heat exchanger, as affected by the various geometrical 
and fluid property factors. Figure 4-7 is for spherical tank with the liquid (Freon-11) 
level at 48.3 cm (19 in.). The tank was locked up and allowed to stratify before the 
coolant was turned on at 46,800 sec. Both the fluid and tank wall temperatures re- 
sponded very quickly, although those nearest the bottom and the top of the tank have the 
fastest response . At the end of the lockup period, there was approximately 1. 7°K 
(3 °R) stratification in the liquid, and an additional 1.7°K (3°R) in the ullage. At the 
end of the cooldown period the spread had increased to 3.6 and 5.6°K (6. 5 and 10°R) 
respectively. The mean cooling coefficient between liquid and wall was approximately 
15% below the theoretical value for the mean Rayleigh Number of 3.15 x 10 which 
would result in a slower pressure response than that predicted by the model. This was 
computed on the basis of experimental mean temperature differences, with boundary 
layer run length as the characteristic dimension. 

Figure 4-8 is for the sphere with a 50% ullage. The same general response is noted; 
i.e. , the temperatures near the top and bottom drop off faster initially than those in 
the center of the tank. The experimental cooling coefficient is approximately 65% lower 
than the theoretical value based on boundary layer run length. However, if liquid, 
height is used in the Rayleigh number, the discrepancy is reduced to 35%. 

A comparison between Figs. 4-9 and 4-7 show the influence of the polyethelene 
spheres which were put into the tank to retard convection. However, the reverse 
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Table 4-5 
TEST SUMMARY 


Test No. 

Apparatus 

Fluid 

% Ullage 

Porous Matrix 

Initial Temperature 
*K <*R) 

Coolant Temperature 
•K (“R) 

Heat Flux 
W/m2 (Btu/Hr) 

A -IGF 

Sphere 

F-ll 

5 

None 

288 

(519) 

278 


5.4 

(1.7) 

A-l 

Sphere 

F-ll 

5 

None 

288 

(519) 

276 


5.4 

(1.7) 

A-2 

Sphere 

F-ll 

15 

None 

292 

(526) 

278 

(499) 

5.4 

(1.7) 

A-3 

Sphere 

F-ll 

50 

None 

293 

(528) 

277 

(498) 

5.4 

(1.7) 

A-4 

Sphere 

F-ll 

5 

3.8 cm (1.5) 

293 

(528) 

277 

(498) 

3.2 

(1.0) 

A-5 

Sphere 

F-ll 

15 

3.8 cm (1.5) 

283 


282 

(507) 

6.4 

(2.0) 

A-6 

Sphere 

F-ll 

50 

3.8 cm (1.5) 

284 

(511) 

279 

(503) 

5.4 

(1.7) 

B-l 

Cylinder 

H 2° 

10 (Sphere) 

None 

277 

(498) 

276 

(495) 

178 

(56) 

B-2 

Cylinder 

h 2 o 

15 

None 

277 

(499) 

276 

(497) 



B-3 

Cylinder 

h 2 o 

5 

None 

278 

(501) 

275 

(495) 

178 

(56) 

B-4 

Cylinder 

V 

50 

None 

* 






B-5 

Cylinder 

F-ll 

10 (Sphere) 

None 

279 

(503) 

276 

(497) 

165 

(54) 

B-6 

Cylinder 

F-ll 

5 

None 

- 

- 


_ 



B-7 

Cylinder 

F-ll 

15 

None 

278 

(501) 

274 

(494) 

165 

(54) 

B-3 

Cylinder 

F-ll 

50 

None 

277 

(499) 

274 

(493) 

86 

(27) 

B-9 

Cylinder 

F-ll 

10 (Sphere) 

1.9 cm (0.75) 

279 

(503) 

276 

(498) 

174 

(55) 

B-10 

Cylinder 

F-ll 

5 

1.9 cm (0.75) 

277 

(499) 

274 

(494) 

(156) 

(49) 

B-ll 

Cylinder 

F-ll 

15 

1.9 cm (0.75) 

280 

(504) 

275 

(495) 

344 

(109) 

B-12 

Cylinder 

F-ll 

50 

1.9 cm (0.75) 

277 

(499) 

275 

(494) 

344 

(109) 

B-13 

Cylinder 

« 2 ° 

10 (Sphere) 

3.8 cm (1.5) 

279 

(502) 

275 

(494) 

344 

(109) 

B-14 

Cylinder 

V 

5 

3.8 cm (1.5) 

277 

(499) 

276 

(497) 

344 

(109) 

B-15 

Cylinder 

H,o 

15 

3.8 cm (1.5) 

275 

(494) 

275 

(494) 

330 

(105) 

B-16 

Cylinder 

F-ll 

10 (Sphere) 

3.8 cm (1.5) 

275 

(496) 

275 

(496) 

330 

(106) 

B-17 

Cylinder 

F-ll 

5 

3.8 cm (1.5) 

280 

(504) 

278 

(501) 

340 

(108) 

B-19 

Cylinder 

F-ll 

15 

3.8 cm (1.5) 

280 

(504) 

277 

(499) 

340 

(108) 

B-19 

Cylinder 

F-ll 

50 

3.8 cm (1.5) 

284 

(511) 

277 

(499) 

340 

(108) 

B-20 

Cylinder 

F-ll 

15 

1. 9 cm (0.75) 


“ 

— 


- 


B-21 

Cylinder 

F-ll 

10 (Sphere) 

1.9 cm (0. 75) 

285 

(513) 

278 

(501) 

338 

(107) 

B-22 

Cylinder 

F-ll 

15 

1.9 cm (0. 75) 

284 

(511) 

278 

(601) 

338 

(107) 

B-23 

Cylinder 

Oil 

10 (Sphere) 

1,9 cm (D.76) 

279 

(502) 

278 

(500) 

338 

(107) 

B-24 

Cylinder 

Oil 

5 

1.9 cm (0.75) 

281 

(506) 

278 

(500) 

340 

(108) 

B-25 

Cylinder 

Oil 

15 

1.9 cm (0.75) 

283 

(510) 

278 

(501) 

338 

(107) 

B-26 

Cylinder 

Oil 

10 (Sphere) 

None 

283 

(510) 

278 

(501) 

338 

(107) 

B-27 

Cylinder 

Oil 

5 

None 

283 

(509) 

278 

(500) 

338 

(107) 

B-28 

Cylinder 

Oil 

15 

None 

282 

(508) 

277 

(499) 

338 

(107) 

B-29 

Cylinder 

F-ll 

10 (Sphere) 

None 

284 

(511) 

277 

(498) 

336 

(107) 

B-30 

Cylinder 

F-ll 

5 

None 

284 

(510) 

276 

(496) 

338 

(107) 

B-31 

Cylinder 

F-U 

10 

None 

284 

(511) 

277 

(498) 

338 

(107) 

C-l 

Sphere 

ln 2 

15 

None 

79.9 

(143) 

65.5 

(HO) 

30.6 

(9.7) 

C-2 

Sphere 

LN 2 

50 

None 

78.8 

(142) 

65.0 

(117) 

30.6 

(9.7) 

C-la 

Sphere 

x^ 2 

15 

None 

76.2 

(141) 

65.5 

(118) 

30.6 

(9.7) 

C-3 

Sphere 

^2 

15 

None 

19.8 

(35.6) 

17.7 

(31.9) 

21.4 

(6.8) 

C-3a 

Sphere 

^2 

15 

None 

20.6 

(37.0) 

18.8 

(33.9) 

21.4 

(6.8) 

C-4 

Sphere 

lh 2 

50 

None 

20.2 

(36.4) 

18.0 

(32.4) 

21.4 

(6.8) 
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appears to be the case. The cooling rates and film coefficients were actually higher 
with the porous media. The theoretical porosity of this media is 0.476 and the pore 
size is the same order as the boundary layer thickness. Consequently, the media may 
have increased the turbulence rather than inhibit growth of the boundary layer. The 
average film coefficient for this run was 15% higher than the thoretical value, computed 
without regard to the media. 

Figure 4-10 is for water in the cylindrical test apparatus. In Test B-3, the mean 
coolant temperature was approximately 274 "K {493° R) until 42,000 sec reaching 282 °K 
(507°F) at 52,000 sec when the test was discontinued. Figure 4-11 shows the wall 
temperatures reversing their trend at approximately 42,000, and leading the reverse 
in the liquid temperatures, as the mode goes from cooling to heating. The cooling 
rate decreased from 17.3 to 10.3 W (59 to 35 Btu/hr) before the coolant temperature 
excursion started a heating trend. The average Rayleigh Number for this test was 

g 

computed to be 2. 5 x 10 which is in the laminar regime. The average cooling coef- 
ficient was 28% higher than the theoretical laminar value. 

Figure 4-12 is for essentially the same conditions as for Fig. 4-10 except that the test 
fluid is oil. The Rayleigh Number for this run was approximately two orders of magni- 
tude lower than for Test B-3. However, it is obvious that the heat exchanger is still 
effective . 

A comparison between Tests B-5 and B-9 (Figs. 4-13 and 4-14) indicate that the 3/4-in. 
media apparently did not reduce the convection as planned. Indeed, the cooling rate for 
Test B-9 was between 46 and 34 W (160 and 116 Btu/hr) with a coolant temperature of 
approximately 273°K (507°R). Pressure was decreasing at approximately 0.35 psi per 
hour. For Test B-5 the cooling rate was between 28 and 11 W (96 and 37 Btu/hr) and 
the rate of pressure change was 0.18 psi per hour. The average coolant temperature 
was 276°K (496°R) for this test. Test B-16 (Fig. 4-15) used the l-l/2-in. media and 
it confirms the fact that low gravity convection w T as not simulated. The cooling rate 
varied from 25 to 6.1W (86 Btu/hr to 21 Btu/hr) with the same coolant temperature 
as for Test B-9. The lower rates are due to the fact that the tank contents are approxi- 
mately 6°K (10°R) lower than for Test B-9. 
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Further confirmation of the ineffectiveness of the media is reflected in Figs. 4-16 and 
4-17. The difference between them is in the ullage configuration which has an insig- 
nificant effect on the energy distribution. Moreover, when compared to Fig. 4-12, one 
notes similar temperature response with and without the media. 

Figures 4-18 and 4-19 primarily illustrate the influence of ullage volume on the per- 
formance of the tank wall heat exchanger. The heat flux for Test B-7 is twice that for 
B-8, but the rate of temperature rise is three times as high and the rates of pressure 
rise are nearly identical. This indicates that the vaporization rates at the interface 
are higher for the large ullage during self-pressurization. However, once the heat 
exchanger was turned on, both the temperature and pressure responses were nearly 
the same. Figure 4-20 is directly comparable to 4-19 with the addition of the 3/4-in. 
media. Once again it reveals no marked effect on performance. 

Figures 4-21 and 4-22 show temperature histories for liquid hydrogen with 15 and 50% 
ullage , respectively. It is apparent that the heat exchanger was not keeping up with the 
incoming heat. The heat rate into the hydrogen tank was determined to be 2.2 W/m 
(6.8 Btu/hr-ft 2 ), which was too high for this heat exchanger design. The maximum heat 
rate that a tank wall heat exchanger can handle is given by: 


'max 


R 


cw 


< t b- t c> 


where ( T - T ) is the mean difference between the bulk fluid and the coolant. 

B C 

For the hydrogen runs, this is the temperature drop across the expansion valve which 
was 4.4°K (8°R) for these tests. The mean value of the tank wall resistance was de- 
termined to be 0.57 °K-m 2 /W (0.4 °R-ft 2 -hr/Btu). Taking into account the reduction 
in thermal conductivity of the 2021 aluminum, the corresponding resistance at liquid 
hydrogen temperatures is 4.85 °K - M 2 /W (3.4 °R-ft -hr/Btu). 
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Consequently, the absolute maximum heat flux to the hydrogen that could be accom- 
modated is: 

q = 4.44/4.85 = 0. 9 W/m 2 ( 0. 28 Btu/hr-ft 2 ) 


which is less than half the experimental value. The pressure histories for these hydro- 
gen tests are shown in Fig. 4-23. 

4.3.2 Comparison With Analytical Convection Model 


In the development of the analytical model, it was assumed that the film coefficients 
between the liquid and the tank wall were given by the conventional correlations for 
vertical flat plate. Example: 



) 


n 


where: 

C 

n 


0.1 

1/3 


| for R a > 10 8 


C 

n 


0.54 

1/4 


for R a < 10® 


In addition, it was assumed that the fluid stratifies when heated, such that the surface 
temperature rise (which establishes the pressure) is given by: 

T g - To 2 (Xc/ r )<£H 
q w (Xc/ R ) = F Pr (A/Xo) 


The temperature at any location X is given by: 
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.2 UU h6 U3 £2 

DYT'EC TIME - 1000 Seconds 

Fig. 4-7 Fluid Temperature Histories for Test A-l, Test Tank-Sphere, Test Fluid- F-ll, Ullage 
Volume 5% 
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DYMEC TIME - 1000 Seconds 

Fig. 4-8 Fluid Temperature Histories for Test A-3, Test Tank-Sphere, Test Fluid- 
F-ll, Ullage Volume 50% 
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Dymec Time - 10 Seconds 


Fig. 4 


9 Fluid Temperature Histories for Test A-4 
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Fig, 4-10 Liquid Temperature Histories, Test B-3 
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Fig. 4-11 Wall Temperature Histories, Test B-3 
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Dymec Time - 10^ Seconds 


Fig. 4-13 Liquid Temperature Histories, Test B-5; Fluid-F-11; 10% Spherical Ullage; No porous Media 
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Fig. 4-14 Liquid Temperature Histories, Test B-9; Fluid-F-11; 10% Spherical Ullage; 3/4-In. Media 
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Fig. 4-15 Liquid Temperature Histories, Test B-16; Fluid F-ll; 10% Spherical Ullage; l-l/2-In. Porous 
Media 
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TIKE FROM START OF TEST- Second 


Fig. 4-16 Liquid Temperature Histories for Test 
B-23, Test Fluid- oil, Test Tank- 
Cylinder, Ullage Volume 10% Spherical 
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Fig. 4-19 Liquid Temperature Histories for Test B-7, Test Fluid-F-11, list Tank - Cylinder 
Ullage Volume 15%, No Internal Porous Media 
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Fig. 4-20 Liquid Temperature Histories for Test 
B-22, Test Fluid - F-ll, Test Tank - 
Cylinder, Ullage Volume 10% Spherical 
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CYI-IEC TIME - 1000 Seconds 


Fig. 4-21 Temperature Histories for Liquid Hydrogen Test C-3; Ullage Volume 15% 
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DYMEC TIME - 1000 Seconds 

Fig. 4-22 Temperature Histories for Liquid Hydrogen Test C-4; Ullage Volume 50% 






TIKE FROM START - IOOO Seconds 


Fig. 4-23 Pressure Histories for Liquid Hydrogen Tests C-3 and C-4 
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Assuming the characteristics shape of the stratified temperature profile to be time in- 
variant, F is given by: 


f T - T. 

F = t ~ _ ~ T '~ d < Z / A > 

J s i 

o 

where Z is the distance below the interface and A is the total depth of the tempera- 
ture stratification. Figures 4-24 and 4-25 show typical temperature profiles at various 
times from the start of heating for Tests B-7 and B-ll. In each run the interface is at 
51.8 cm (20.4 in.). Test B-7 has only freon in the tank, while B-ll has Freon with the 
3/4-in. polethelene media. The respective heat fluxes are 136 and 334W/m , (43 and 

o 

106 Btu/hr-ft ). The stratified layer depth was obtained and the nondimens ional pro- 
files were plotted as shown in Figs. 4-26 and 4-27. These confirm the time invariancy 
of the profile shape and result in energy integrals of 0.42 and 0.37 for Tests B-7 and 
B-ll, respectively. The dashed lines shown in Fig. 4-26 are obtained with the above 
equation for ( T - T B ) using the experimental value F =0.42. 

Figure 4-28 shows the temperatures for Test B-26. The test fluid is oil. Also, the 
tank is filled all the way to the top; the ullage is spherical; and the 3/4-in. porous media 
is used. Even discounting the porous media, the effective Rayleigh Number is 6 orders 
of magnitude below that for Run B-7. Neverthless, the value of F is 0.4 as shown on 
Fig. 4-29. Figure 4-30 shows the nondime ns ional profile for Test B-l (water) for which 
F is approximately 0.4 also. 

Figure 4-31 shows some temperature profiles for Test A-2 during warm-up period. 
These are nondimensionalized on Fig. 4-32. Again, the time invariancy is confirmed. 
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Fig. 4-24 Temperature Stratification Profiles, Test B-7 
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Fig. 4-25 Temperature Stratification Profiles, Test B-ll 
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Fig. 4-26 Nondimensional Temperature Profiles, Test B-7 
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Fig. 4-27 Nondimensional 
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Fig. 4-28 Temperature Profiles for Test B-26, Test Fluid-Oil, Ullage Volume 5%, 
3/4-in. Porous Media 
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Fig. 4-30 Nondime nsional Temperature Stratification Profile, Test B-l 





























Fig. 4-32 Nondimensional Temperature Stratification Profile, Test A-2 



The value of F , determined by graphical integration of Fig. 4-32, is 0.37. The liquid 
level for this run was at approximately 17 in. which corresponds to the maximum liquid 
temperature. The tank was vented just before the beginning of the test which caused 
some instability in the vapor temperatures above the interface. Thus, the temperature 
at 50.9 cm (20 in.) started off slightly subcooled relative to the surface liquid. However, 
within 1200 sec, the vapor exceeded the liquid temperature and a stable condition 
existed before the heat exchanger was activated. 

Figure 4-33 compares the results from these tests with the correlation for F which had 
been developed and reported in Ref. 8. The data from this program are in good agree- 
ment with the correlation. 

The other major assumption employed in the model is the relationship between the 
Nusselt and Rayleigh Numbers. The average heat transfer coefficients between the 
liquid and wall were computed from the measured cooling rates and temperatures at 
several times during a test. Nusselt Number was then calculated, using boundary layer 
run length as the characteristic dimension, for comparison with a theoretical value 
determined from the experimental Rayleight Number. The temperature difference used 
to determine both the experimental film coefficient and the Rayleigh Number was the 
difference between the mean liquid and mean wall temperatures, i.e. , 



For a few cases, this was compared with the mean of the differences, i.e. , 

at wl = 4 2 ( t l - t w) l 

which was evaluated graphically. This method was considerably more time-consuming 
and did not improve the data correlation. Therefore, it was discarded in favor of the 
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Figure Comparison Between Energy Integral Data Form Current Program and That of Previous 

Investigators 

Fig. 4-33 Comparison Between Energy Integral Data From Current Program and That of 
Previous Investigators 


difference in the mean temperatures which could be determined rapidly with the aid of 
an electronic calculator. 

The results from this portion of the comparison between model and experiment are 
shown on Fig. 4-34. The data generally support the model, although there is consider- 
able scatter, particularly at the high Nusselt Numbers. However, the Nusselt temper- 
ature differences were rather small (0.2°K — 5°K). The precision of the thermocouples 
was estimated to be approximately 0.3°K (0.5°R). It is believed that much of this 
scatter is due to these small temperature differences. 

Figure 4-34 also shows the ineffectiveness of the polyethelene spheres. If the matrix 
had been effective in simulating low-g connection, the Nusselt Numbers for those tests 
(solid symbols) should be nearly two orders of magnitude lower than for the other 
tests, which is clearly not the case. Additional work in this area should be consider- 
ed wherein smaller apparatus and much smaller pore sizes are used. 

The confirmation of the heating profile and the cooling coefficients lend a certain 
amount of validity to the analytical model. However, the real proof comes not from 
the individual parts, but rather from the integrated model and the extent to which it 
correlates the test data. Figure 4-35 shows a direct comparison between the model 
prediction and the experimental pressure response for Test B-7 which was conducted 
in the cylindrical tank with Freon -11. Figures 4-36 and 4-37 show the temperature' 
histories for the same test. The correlation is very good. On Figure 4-37, the dashed 
line represents the test data corresponding to the upper liquid node (node 6) and the 
bottom liquid node (node 14). At intermediate locations, the temperatures are nearly 
parallel. It can be seen that the major difference between model prediction and test 
data is that the computed temperatures merge into the surface temperature during the 
cooling process, which is not the case with the experimental data. The temperature 
merging occurs in the model when tank pressure drops to local saturation pressure, 
because no allowance is made for the static head on the subsurface liquid. 
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In the development of the analytical model, it was assumed that a spherical tank could be 
represented as a psuedo- cylindrical tank for which the geometry, psuedo heating rates, 
etc. , are computed in the GEOM subroutine based on spherical tank inputs. The validity 
of this can be assessed from Figs. 4-38 and 4-39. These present a direct comparison 
between experimental and pressure responses with liquid hydrogen at 15 and 50% ullage. 
The curves indicated as "model” are based upon the computation with a uniform heat 
flux established from the boiloff tests. However, when the tank wall heat exchanger was 
activated, and the tank was locked up, the vent line became a trapped vapor pocket which 
extended from the top of the tank out to the chamber wall. The effect on ullage volume 
was insignificant. However, this results in a heat leak that is not determined during a 
boiloff test when the line is being cooled by the vented vapor, but which goes into direct 
ullage heating during the tank wall heat exchanger operation. An upper and lower limit 
was estimated to be 0.85 W (2.9 Btu/hr) and 0.22 W (0.75 Btu/hr), respectively, for 
this direct heat source. The curves labeled as "modified model” on Figs. 4-38 and 
4-39 are for the added effect of the direct addition of 0.85 W to the ullage. The basic 
model falls below the experimental pressure resp nses and the upper limit overpredicts 
the test results. Therefore, the deficiency is in the experiment and one can conclude 
that the analytical model developed in this program, for the tank wall heat exchanger, is 
a valid performance evaluation instrument. 
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RAYLEIGH NUMBER 

Fig. 4-34 Comparison Between Experimental Nusselt Numbers and Analytical Model 
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Fig. 4-35 Comparison Between Experimental Pressure History and Analyti- 
cal Model Prediction, Test B-7 
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Fig. 4-36 Comparison Between Experimental and Analytical Saturation 
Temperatures, Test B-7 
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Fig. 4-37 Comparison Between Experimental and Analytical Liquid Temperatures 
Test B-7 ' 
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Fig. 4-39 Comparison Between Experimental Pressure and Analytical Prediction 
For Test C-4; Test Fluid-LH2; Ullage Volume 50% 
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Section 5 

CONCLUDING REMARKS 


This program was conducted to develop an analysis technique for evaluating the 
performance of a tank wall heat exchanger as a means of controlling pressure in a 
cryogenic tank during space flight. As a result of this work, the following conclusions 
can be drawn: 

(1) A tank wall heat exchanger can be used to control tank pressure. 

(2) Some propellant subcooling does occur, resulting in excess loss of propel- 
lant. However, the net effect of this inefficiency is insignificant after 5 or 
6 vent cycles . 

(3) The effectiveness of the tank wail heat exchanger is not diminished in the 
absence of gravity because subsurface boiling tends to mix the propellants. 

(4) A heat exchanger mounted on the external wall is heavier than one mounted 
internally, but is probably preferred from a reliability standpoint. 

(5) Tests conducted in this program have confirmed the validity of the analytical 
model. Therefore, the parametric data that have been developed with the 
model can be used to estimate the pressure control characteristics for the 
range of geometries and environments covered herein. 

(6) A packed sphere porous media does not appear to retard free convection for 
the size of test articles used in this program. Additional investigations for 
different porous media, and with smaller test articles, should be conducted 
to extend the results from this program down to lower Rayleigh Numbers 
(lower gravity simulation). 


PRECEDING PAGE BLANK NOT FILMED 
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Appendix A 
TANK FILL MODEL 


o During the first phase of the fill process liquid which enters the tank is 
saturated and vaporized instantaneously, mixing with, and being raised to 
the temperature of, the bulk of the ullage vapor* 

o The first phase of fill continues until the ullage vapor temperature has 
dropped to the point at which a liquid pool forms at the one end of the 
tank. This point in time is determined by testing the evaporation time 
of an "average" droplet. 

o During the second phase of fill, the liquid forms a pool at one end of the 
tank with incoming liquid going directly into this pool and the ullage con- 
sisting only of vapor (no liquid droplets). 

^ * 

Phase I Fill (T > T ) 
u. s 

With the above assumptions, the general receiver tank equations for phase I are: 



(A-l) 


(A-2) 



(A -3) 


M* T* 
u u 

P* 


1 


(A-4) 


A-l 
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Equations (2) and (4) give 


dT 
u 

dt* 


1 dp* 

M * dt* 
u 


p* 

M *2 


dM * 
u 

dt* 


r 



(A -5) 


Substituting (5) into (3) and using (4) yields 



Solving Eq. (1) for T* , 

s 


= 1 - B t / RTi p *\ 

T* C t A £ n(P.v P ; 
1 ' 


(A -7) 


The ullage vapor temperature is obtained from Eq. (4) as 


u 


p* 

M* 

u 


(A -8) 


The ullage mass is obtained from 
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or 



m * ) dt* 
v ' 


(A -9) 


Substituting (9) into (8) yields 


rp* 

u 



*p * > T + 

u — s 


(A -10) 


The ullage wall temperature is obtained from a heat balance on the ullage wall 


dT* 


uw 

dt* 


i °i 


M* 

uw 


C 

uw 


Q* 

quw 


h* A* (T* 
uw uw uw 


T *) 

u ’ 


(A -11) 


Equations (6) and (11) must be solved simultaneously, with the aid of Equations (7) 
and (10). 
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Phase II Fill <T *«T *) 

x u s 

The simplified model assumed is not valid when the gas temperature nears the saturation 
temperature, since there is not sufficient heat transfer to vaporize all the liquid as it 
enters the tank. For this reason, Eq. (6) is integrated only until the ullage vapor 
temperature is approximately equal to the saturation temperature. It is at this point 
at which the Phase II analysis is applied. 


For this analysis it is assumed that: 

o Gas phase only exists in the ullage (M|y=0) 

o Liquid entering the tank is subcooled having a dimensionless temperature 


T+ and X=0. 

o No boiling at tank -wall bottomed-liquid interface (rhj sw =0) 
The general receiver tank equations then reduce to 

3 

RT. 


p* = 


on s tnen reduce to 

-tM] 


(A —12) 


dM* 

Jj 

dt* 


dM* 
u 

dt* 





= *r* - T [Vs ( T u - T l) - h sl ( T s - t b)] 
= ^r[ h us (k - T t)- h s L ( T s - t b] - 


(A -13) 


(A -14) 


M bS- = “f ( T L - T b)- [ h us ( T ; - T s)- ( T S - T b)} + A S h *SL 

(t* - T *\ + A* q* 
y L S a B ) A WL q BL 


where qg L is the wall heat added to the liquid. 



(A - 16 ) 
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M* T* 

M* + — t — = 1 
B p* 


(A -17) 


Equations (14) and (17) yield 


dT* * dM* (l-M*A , * 

u_ -p* B \ B / dp* 

dt* M* dt* M* dt* 

u u 


Substituting (18) into (16) and using (17) yields 


dp* _ 1 I u u ^ , rp* 

dt* ~(1-M*) ((l-M^) dt* u dt* 


dM* 

M*2 dt* 
u 


(A-18) 


(t* - 

V u 

T *\ 

— -A* 

h* | 

f T* - T* ^ 

I" h CT 

It* - t*\] 

A 

s 

1 us 1 

l u sj 

1 SL 

\ 8 B )J 


pv 


fh* A*/T*-T*\- h* A* fa* -T*\l( 
us s \ u s) wu wu\ wu u m 


(A-19) 


The interface temperature is obtained from (12): 

'p* — ' 


- cfr H 


'RT 


The masses M* and M* are obtained from integration of equations (13) and (14); 

t* 

dM* \ 

dt* 


i 


and 


t* 


m ; * /(£?}*• 
t* X ' 
l 


(A -20) 


(A -21) 


(A -22) 


Where t? is the dimensionless time at the beginning of Phase n fill, Mu* is the ullage 
mass at t*, and M* = 0 at the t* . 

1 D 1 
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The ullage vapor temperature is obtained from (17): 


= ( 1 - m b) 

The bottomed liquid temperature is derived from integration of Eq. (15): 

t* 



where we assume that T^ = T* . 


(A -23) 


(A -24) 


(A -25) 


The ullage wall temperature is obtained from a heat balance on that portion of the wall 
exposed to the ullage vapor: 


where 


dT* 

uw 

dt* 


M* C 
uw uw 


[■ 


a* - h* A* 
niw uw uw 


- 

/T* - T*\ 
\ uw u) 


q* = 
^uw 


uw 


h 1 L 2 T 
o o o 


* Saw ^o 


w 


uw 


w 


i* =■ 

*BL 


*BL 

h 1 L 2 T 
o o o 


*BL 


= q 


BL 


o A„ 
w w 


= q. 


(a - a ) 

V w uw / 


w 


w 


(A-26) 


Criterion for Transition From Phase I to Phase II 

— -3 

Assume mean drop diameter D*^ - 1. 73 x 10 ft. 


From LMSC-A847685, P.A-27, droplet will evaporate before traversing ullage if 

2 / \ / \ 


, K u<V T s>V R I 
aT \ D D 


-M> 10" 3 


1 


where Rj is radius of inlet. 


= U _ x )+ X 1 = for X = 0, so 

x |_ p ? v 7 J p u 


* rsr t 
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Therefore, transition from Phase I to Phase II occurs when 

Xm f Dj^xKf 3 


or 


For Rj in feet. 


( T u - T s )< 
( T u- T s)< 


K R/ 
u I 




K R/ 
u I 


x 10 


-9 


(A- 


Re suits 


The model was used to calculate pressure histories for a tank that has been depleted 
and is being refilled in orbit. It was assumed that the tank was not evacuated prior to 
the fill and that both the tank and its contents were at 244°K (440°R)« Starting pressure 
was assumed to be 1 atmosphere. Figure A-l is a typical result. It was found that 
there is a very small drop in pressure (» 0. 1 psi) when the cold fluid first enters the 
tank. Thereafter, the pressure starts to climb rapidly until the ullage and saturation 
pressure approach equilibrium, near the end of Phase 1. Apparently, the combination 
of mass addition and heating from the tank walls more than offsets the cooling of the 
original gas. After the walls have been cooled to near saturation the pressure begins 
to decline slowly. Flow rates were varied from 0. 65 percent to 65 percent of tank 
capacity per hour with only a small effect on the peak pressure. However, the time at 
which peak pressure occurs is nearly inverse to the fill rate. On the basis of this 
analysis, it appears that if a tank is to be refilled in space, it should first be evacuated 
and cooled with the vent open. 
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Fig* A-l Typical Pressure History During an In Space Refilling Operation 
With LH2; Tank Diameter = 1.2m (4 Ft), Fill Rate = 6.5 Percent 
Tank Capacity per Hour 
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DEFINITION OF SYMBOLS 


2 2 

Cross-sectional area for flow, ft (m ) 

2 2 

Interfacial surface area between bottomed liquid and tank' ullage, ft (m ) 

Interfacia I surface area between dispersed-liquid system and ullage-vapor 
system, ft* ( m 2) 

Specific heat at constant pressure, Btu/lbm-°F (J/kg-°K) 

Specific heat at constant volume, Btu/lbm-°F (J/kg-°K) 

Droplet diameter, ft (m) 

Conversion ratio, 32.2 lbm-ft/lbf-s^ 

2 2 

Acceleration of gravity, ft/s (m/s z ) 

2 2 

Heat-transfer coefficient based on surface area, Btu/hr-ft -°F (W/m °K) 
Mechanical equivalent of heat, J = 778 ft Ibf/Btu 
Thermal conductivity, Btu/hr-ft -°F (W/m °K) 

Mass, Ibm (kg) 

Mass flowrate, Ibm/hr (kg/hr) 

Mass flowrate of fluid entering tank through fill line, Ibm/hr (kg/hr) 

Mass rate of vapor addition to tank ullage due to vaporization at the free 
surface minus the mass rate of condensation of incoming vapor which con- 
denses directly into the bottomed liquid, Ibm/hr (kg/hr) 

Mass flow rate of vapor through tank vent, Ibm/hr (kg/hr) 

Pressure, lb/ft^ or psi (N|/cm^) 

Volumertric flow rate, ft^/hr (m 3/5) 

Heat flux density, Btu/hr-ft^ (W/m^) 

2 2 

Heat flux density into bottomed liquid from tank wall, Btu/hr-ft (W/m ) 

0 o 

Heat fl ux density into ullage vapor from tank wall. Btu/hr-ft^ (W/rn ) 

Gas constant, ft-lbf/lbm-°R (J/°K mol) 

2 2 

Surface area, ft (m ) 

Temperature, °F (°K) 
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DEFINITION OF SYMBOLS (Cont.) 


Time, hr 

Volume, ft^ (m^) 

Specific volume, ft^/lbm (m^/kg) 

3 3 

Specific volume change during evaporation, frylbm ( m /kg) 
Vapor compressibility factor 

Greek Letters 

Constant, ratio of specific heats 

f < : t 

Void fraction or porosity 
Constant 

Dimensionless heat of vaporization 
Latent heat of vaporization, Btu/lbm ( j/kg) 

Density, Ibm/ft^ (kg/m^) 

<x = 0!l714x l0‘ 8 Btu/hr-ft 2 -°R 4 ' 

= 5.6697 x 1 O' 8 W/m 2 °K 4 , 

- ! ’ - "■>' 

A function 

Superscripts 

Refers to nondime ns iona I quantity 
Refers to time average or mean value 
Refers to critical value 

•i v 1 > ’ . 

Subscripts 

Bottomed-liquid condition . .„■* 

Bubble property or bulk -flow condition 
Core condition 

Droplet or dispersed-liquid condition 
Drop condition at tank inlet 
Exit condition 
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DEFINITION OF SYMBOLS (Cont.) 


f 

9 

St 

LI 

o 

s 

T 

u 

V 

v 

w 


Fill condition 

Gas or gaseous-phase condition 
Inner diameter or surface 
Saturated liquid or liquid phase 
Inlet-liquid condition 
Initial condition or outer surface 
Saturation condition * 

Tank Condition 
Ullage vapor condition 
Vapor condition 
Vent condition 
Wall condition 
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Dimensionless Groups 



L q “ Characteristic tank dimension (diameter for spherical tank)* 
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Appendix B 

HEAT EXCHANGER PRESSURE DROP 


The liquid extracted from the tank is expanded to a low pressure and temperature and 
passed through the heat exchanger to absorb enough energy to vaporize the two-phase 
mixture. This vented vapor is then vented into the vacuum of space. 


From a momentum balance, as the fluid flows through the heat exchanger, the pres- 
sure gradient at any point in the tube is given by: 



where the subscripts, T , f , a , g refer to the total, frictional, acceleration, and 
gravitational gradients, respectively. 


Dr. Owens* has shown that the frictional gradient with two-phase flow is related to 
that with liquid only flowing in the tube, by the following expression: 



where x is the mass fraction of vapor, and subscripts g , i pertain respectively to 
vapor and liquid properties. The frictional drop with liquid is given by: 



n g 2 vf 

2Dg c 


*W. L. Owens, "Two-Phase Pressure Gradients" International Developments in Heat 
Transfer, Part II, pp. 363-368, January 1962 
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The acceleration pressure gradient is given by: 


/dP \ = dv 

\ d i/ a S c d l 


where: 


V : 

II 

-x) 

+ X V 

g 


dv 

. av 9 

v i 


dx 3v 


3v i 

di 

+ dx 

di ~W g 



d v 


= (1 “ x) 


i 


(if). 


+ (V 


g 



The gradient due to the gravity vector is given by: 



g/g. 


sin 9 
v 


= g/g. 


sin 9 

(1 “ X) + X V g 


Substituting these various expressions into the first equation, we get for the total 
pressure gradient at any point 





dx 

V i di + 


g/g sin e 


A /I 


1 4 - 


gM 

e c \ 


(1 -x) 


dv J? 

dP 


d v g | 
dP j 


d v» 

For cryogens, the term (1 - x) - r g is insignificant compared to the vapor term 
d v 

£ 

x -^p- 1 . Therefore, the preceding expression can be somewhat simplified to give: 



Numerator 


1 + x 


r 2 d V rr 

9 E. 

g. dP 
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The total pressure drop in the heat exchanger tube must be obtained by integration, 

i* 6* f 



This equation is not amenable to direct solution except by the tedium of trial and error 
with step wise integration. However, we note that in the denominator, the quantity, 
r 2 d v 2 

— ■■ A i s the (ratio) of actual to sonic velocity if the vent fluid is all vaporized. 

d v 

It must be less than 1 in the heat exchanger tube. Since is a function of 

pressure and the mass flux is determined by the vent rate and tube size, we can use 
the limit condition to estimate the minimum tube size that can be used. Fig. B-l 
shows the maximum mass flux as a function of pressure in the tube. 


The mass flux in the heat exchanger is related to the cooling rate and the tube diameter, 

i. e. , 


ft 

v 

^tube 


4 q 

c 

Xir d 


2 

t 


where: 

& = vent flow rate 

v 

A = tank surface area 
s 

dj. = tube diameter 

q - average cooling flux at the tank wall 

G 
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Also, the cooling flux is related to the average wall resistance, R , by: 

' , * k .} . GW 


q 


c 


" » > < T L ' T c> 
R 

cw 


Combining the two preceding equations and replacing tube diameter d^ by tab diam- 
eter 2a , we get: 


G 



(1 - n ) (T t - t ) 

AJ 

^ R cw < a > 2 


This equation can be rearranged to describe the minimum tube size in terms of tank 

size, R and the maximum mass flux, 

* cw 


a 

mm 



(1 -r» (T b - T c > 


R. 


cw 


7I\G 


max 


1/2 


The effectiveness has been shown to be related to the value of R cw . The expansion 

temperature drop, (T-^.- T ) is related to. the heat exchanger pressure, for a given 

. — 

tank operating pressure, as is the maximum mass flux. Therefore, for a given 
pressure at the tube inlet, the last equation can be used to estimate the minimum tube 
size for any tank and design value for R . 

CW 
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APPENDIX C 


List of Symbols 

Heat exchanger attachment radius 
A - Area 

A_ - Flat interface area 
F 

B 

B - Bond Number 
o 

C - Sensible heat capacity 

Cp - Specific heat at constant pressure 

Cy - Specific heat at constant volume 

D - Tank diameter 

f - Friction factor inside heat exchanger tube 

F - Energy integral defined by temperature profits 

G - Coolant mass flux through heat exchanger 

g - Gravitational acceleration 

g - Standard gravitational acceleration 
o 

g - Universal gravitational constant 

G - Grashof Number 

R 

G* - Modified grashof number based on heat flux 
R 

h - Heat transfer film coefficient 

k - Thermal conductivity 

K - Slope of P-T saturation line 
s 

L - Boundary layer length 

M - Fluid quantity 

- Mass flow rate from boundary layer 

C-l 


M (Ft ) 

M 2 (Ft 2 ) 

M 2 (Ft 2 ) 

dimensionless 

cal/gr ( Btu /lb) 
cal/gr-ok ( B ^ u /lb-°R) 

cal/gr-ok ( B "^ u /lb-°R) 

M (Ft) 

dimensionless 
Kg/sec-M 2 (lb/sec-Ft 2 ) 
M/sec 2 (Ft/sec 2 ) 
M/sec 2 (Ft/sec 2 ) 

dimensionless 

dimensionless 

W / M -° K (HrFt° F ) 

N/M 2 -°K (PSi/°R) 

M (Ft) 

Kg (ibm) 

Kg/hr . ( lbm/hr . ) 
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M 

e 

- 

Mass flow rate into boundary layer 

- Kg/hr . ( lbm/ hr . ) 

N 

- 

Number of heat exchanger attachment points 


Nu 

_ 

Nusselt number 

- dimensionless 

P 

- 

Pressure 

- N/M 2 (lb /Ft 2 ) 

P 

_ 

Prandtl Number 

- dimensionless 

r 



/Btu \ 

Sf 

- 

Heat flux at the tank wall 

- ^(hr-Ft 2 ) 

Q 

- 

Heat rate 

- watts (Btu/Hr) 

% 

- 

Heat rate into each heat exchanger attachment 

- watts (Btu/Hr) 

R 

a 

- 

Rayleigh Number 

- dimensionless 

R* 

a 

- 

Modified Rayleigh Number 

- dimensionless 

R 

- 

Tank radius 

- M (Ft) 

R 

_ 

Thermal resistance between the tank wall 

- M 2 -°K /hr-Ft 2 -°R) 

ew 


and coolant stream 

' Btu ’ 

r f 

- 

Radius at flat interface 

- M (Ft) 

S 

- 

Spacing between heat exchanger attachment 
points 

- M (Ft) 

T 

- 

Temperature 

- °K(°R) 

t 

- 

Tank wall thickness 

- M (Ft) 

t 

- 

Time 

- sec 

U-. 

_ 

Characteristic fully developed, turbulent 

- M/sec (Ft/sec) 

i 


boundary layer velocity 

U 

s 

- 

Average, steady state boundary layer velocity 

- M/sec (Ft/sec) 

V 

- 

Volume 

- M 3 (Ft 3 ) 

V 

- 

Specific volume 

- M 3 /Kg (Ft 3 /lb) 

Wt 

- 

Heat exchanger tube weight 

- Kg (lbs) 

Ws 

- 

Tank wall membrane weight 

- Kg (lbs) 

X 

- 

Axial distance measured from bottom of tank 

- M (Ft) 

y 

- 

Radial distance measured from tank wall 

- M (Ft) 

z 

- 

Compressibility factor 

- dimensionless 

z 

- 

Elevation of flat liquid interface 

- M (Ft) 

O' 

- 

Number of standard gravities 

- dimensionless 

(3 

- 

Liquid coefficient of thermal expansion 

- °k" 1 (°r' 1 ) 
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6 - Boundary layer thickness 

A - Stratified layer thickness 

O' - Surface tension 

Porosity 

r? - Heat exchanger effectiveness 

y - Kinematic viscosity 

Density 

- Latent heat of vaporization 

ip - Geometry factor used in GEOM subroutine 

- Characteristic time for stratification 


- M (Ft) 

- M (Ft) 

- Kg/m (lbp/Ft) 

- dimensionless 

- Kg (lb m ) 

- cal /Kg ( Btu /lb) 

- dimensionless 

- dimensionless 


Subscripts 

B - Pertains to area or volume subjected to 

bottom heating 

C - Pertains to cylindrical geometry 

cm - pertains to cooling at tank wall 

S ~ gas 

L - Liquid 

S - Saturation conditions at liquid surface 

s - Steady state 

s,t - Steady state, turbulent 

s,L - Steady state, evaluated at axial distance L 

T - Total 

TURB - Turbulent 

U - Ullage 

W,t - Pertains to conditions at the tank wall, 
turbulent boundary layer 
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Rocketdyne Division 

Worth American Rockwell, Inc. 

6633 Canoga Avenue 

Canoga Park, California 91304 

Attn: Library, Department 59&-306 

Space & Information Systems Division 
Worth American Rockwell 
12214 Lakewood Blvd 
Downey, California 
Attn: Library 

Northrop Space Laboratories 
3401 West Broadway 
Hawthorne, California 
Attn: Library 

Purdue University 
Lafayette, Indiana 47907 
Attn: Library (Technical) 

Radio Corporation of America 
Astro-Electronics Products 
Princeton, New Jersey 
Attn: Library 

Rocket Research Corporation 
Willow Road at ll 6 th Street 
Redmond, Washington 98052 
Attn : Library 

Stanford Research Institute 
333 Ravenswood Avenue 
Menlo Park, California 9^-025 
Attn: Library 
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R. W. Hallet 
G. W. Burge 
P. Klevatt 


R. A. Herzmark 


Dr. R. J. Thompson 
S. F. Iacobellis 


Dr. William Howard 


Dr. Bruce Reese 


F. McCullough, Jr. 


Dr. Gerald Marksman 
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Copies 


1 


1 


1 


1 


1 


1 


1 


1 


Thiokol Chemical Corporation 
Redstone Division 
Huntsville, Alabama 
Attn: Library 

TRW SysteEis Inc. 

1 Space Park 

Redondo Beach, California 90278 
Attn: Tech. Lib. Doc. Acquisitions 

TRW 

TAPCO Division 
23555 Euclid Avenue 
Cleveland, Ohio 44117 

United Aircraft Corporation 
Corporation Library 
400 Main Street 

East Hartford, Connecticut 06108 
Attn: Library 

United Aircraft Corporation 
Pratt & Whitney Division 
Florida Research & Development Center 
P.0. Box 2691 

West Palm Beach, Florida 33402 
Attn: Library 

United Aircraft Corporation 
United Technology Center 
P.0. Box 358 

Sunnyvale, California 9^038 
Attn: Library 

Vickers Incorporated 
Box 302 

Troy, Michigan 

Vought Astronautics 
Box 5907 
Dallas, Texas 
Attn: Library 


John Goodloe 


D. H. Lee 


P. T. Angell 


Dr. David Rix 
Erie Martin 
Frank Owen 
Wm, E. Taylor 


R. J. Coar 
Dr. Sehmitke 


Dr . David Altman 
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Rockwell International 
Space Division 
12214 Lakewood Boulevard 
Downey, California 90241 

General Dynamic s/Convair 
P. 0. Box 80847 
5001 Kearny Villa Road 
San Diego, California 92138 


Dr . Don Gluck 
Mail Stop SK-06 


Mr. John Stark 
Mail Zone 632-00 


Minnesota Mining & Manufacturing Company 
900 Bush Avenue 
St. Paul, Minnesota 55106 
Attn: Library 
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